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Abstract 

We consider the minimization of J Q |Du| p dx in a perforated domain fL := B E {an) of 

R™. among maps u € W 1,v {Vl E ,W l ) that are incompressible (dctZ?u = 1), invertible, and satisfy 

a Dirichlet boundary condition u = g on <9Q. If the volume enclosed by g(dil) is greater than 

|f2|, any such deformation u is forced to map the small holes B e {sn) onto macroscopically visible 

cavities (which do not disappear as e — > 0). We restrict our attention to the critical exponent 

p = n, where the energy required for cavitation is of the order of X}i=i ^1 log ^1 and the model 

is suited, therefore, for an asymptotic analysis (v\, . . . , i>m denote the volumes of the cavities). 

In the spirit of the analysis of vortices in Ginzburg-Landau theory, we obtain estimates for the 

p 

"renormalized" energy i f Q J^t 1 dx — v i \ log e | , showing its dependence on the size and 
the shape of the cavities, on the initial distance between the cavitation points ai,. . . ,slm, and 
on the distance from these points to the outer boundary dfl. Based on those estimates we 
conclude, for the case of two cavities, that either the cavities prefer to be spherical in shape and 
well separated, or to be very close to each other and appear as a single equivalent round cavity. 
This is in agreement with existing numerical simulations, and is reminiscent of the interaction 
between cavities in the mechanism of ductile fracture by void growth and coalescence. 

1 Introduction 

1.1 Motivation 

In nonlinear elasticity, cavitation is the name given to the sudden formation of cavities in an 
initially perfect material, due to its incompressibility (or near-incompressibility) , in response to 
a sufficiently large and triaxial tension. It plays a central role in the initiation of fracture in 
metals [351 E21 ES EH EE] and in elastomers [291 E3 E21 EH Q"8] (especially in reinforced elastomers 
EHEDQSJEIEI]), via the mechanism of void growth and coalescence. It has important applications 
such as the indirect measurement of mechanical properties [45] or the rubber-toughening of brittle 
polymers |461 [T4"| 1761 [4*8] . Mathematically it constitutes a realistic example of a regular variational 
problem with singular minimizers, and corresponds to the case when the stored-energy function of 
the material is not VF 1,p -quasiconvex [2JEJ|~7], the Jacobian determinant is not weakly continuous 
[7], and important properties such as the invertibility of the deformation may not pass to the weak 
limit |551 Sect. 11]. The problem has been studied by many authors, beginning with Gent-Lindley 
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[30] and Ball @]; see the review papers [29j[HJ[25], the variational models of Miiller-Spector [55] and 
Sivaloganathan-Spector [70], and the recent works [38\ 149] for further motivation and references. 
The standard model in the variational approach to cavitation considers functionals of the form 

\Du\ p dx, (1.1) 

n 

where the deformation u : $7 C W 1 —> M. n is constrained to be incompressible (i.e. detDu = 1) and 
globally invertible, and either a Dirichlet condition u = g or a force boundary condition is applied. 
Unless the boundary condition is exactly compatible with the volume, cavities have to be formed. 
If p < n this can happen while still keeping a finite energy. A typical deformation creating a cavity 
of volume oj n A n at the origin (oj n being the volume of the unit ball in W 1 ) is given by 

x 



u(x) = ^4™ + |x|™ — . (1.2) 



We can easily compute that 



. l2 (n-l)A 2 , . 

Du 2 ~ K > . 1.3 
x=o xr 



In that situation the origin is called a cavitation point, which belongs to the domain space, and 
its image by u is the cavity, belonging to the target space. Contrarily to the usual, we study the 
critical case p = n where the cavity behaviour (jl.2p just fails to be of finite energy. 

This fact is analogous to what happens for § 1 -valued harmonic maps in dimension 2, which 
were particularly studied in the context of the Ginzburg-Landau model, see Bethuel-Brezis-Helein 
HUj. For Si-valued maps u from l!cl 2 , the topological degree of u around a point a is defined 
by the following integer 

, i r du 

a = — / — x u. 



27r JdB{a,r) 9t 

Points around which this is not zero are called vortices. Typical vortices of degree al look like 
u = e tde (in polar coordinates). If d ^ again |Z)u| 2 just fails to be integrable since for the typical 
vortex |-Du| 2 ~ x=0 Mj^ 5 just as above (|1.3|) . up to a constant factor. So there is an analogy in 
that sense between maps from Q to C which are constrained to satisfy |u| = 1, and maps from 
$7 to M. 2 which satisfy the incompressibility constraint det-Du = 1. We see that in this analogy 
(in dimension 2) the volume of the cavity divided by it plays the role of the absolute value of the 
degree for Si-valued maps. In this correspondence two important differences appear: the degree is 
quantized while the cavity volume is not; on the other hand the degree has a sign, which can lead 
to "cancellations" between vortices, while the cavity volume is always positive. 

In the context of Si-valued maps, two possible ways of giving a meaning to |-Du| 2 are the 
following. The first is to relax the constraint |u| = 1 and replace it by a penalization, and study 
instead 

/ |Du| 2 + l(l-|u| 2 ) 2 (1.4) 
Jn e 2 

in the limit e — > 0; this is the Ginzburg-Landau approximation. The second is to study the energy 
with the constraint |u| = 1 but in a punctured domain £l e := 0\ Uj B(a.i,e) where a^'s stand for 
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the vortex locations: 

min / \Du\ 2 (1.5) 

H=iyn s 

again in the limit e — > 0; this can be called the "renormalized energy approach". Both of these 
approaches were followed in [lOj , where it is proven that the Ginzburg-Landau approach essentially 
reduces to the renormalized energy approach. More specifically, when there are vortices at aj, \Du\ 
will behave like |dj|/|x — a,| near each vortex (where d{ is the degree of the vortex) and both energies 
(|1.4p and (|1.5p will blow up like ir d\ log \ as e — > 0. It is shown in [10] that when this divergent 
term is subtracted off (this is the "renormalization" procedure), what remains is a nondivergent 
term depending on the positions of the vortices aj and their degrees di (and the domain), called 
precisely the renormalized energy. That energy is essentially a Coulombian interaction between the 
points aj behaving like charged particles (vortices of same degree repel, those of opposite degrees 
attract) and it can be written down quite explicitly. 

Our goal here is to study cavitation in the same spirit. A first attempt, which would be the 
analogue of (11. 4j) . would be to relax the incompressibility constraint and study for example 

f | DU |* + < 1 - det ° u > 2 . (1.6) 
Jn £ 

We do not however follow this route which seems to present many difficulties (one of them is 
that this energy in two dimensions is scale invariant, and that contrarily to (jl.4p the nonlinearity 
contains as many order of derivatives as the other term), but it remains a seemingly interesting 
open problem, which would have good physical sense. Rather we follow the second approach, i.e. 
that of working in punctured domains while keeping the incompressibility constraint. 

For the sake of generality we consider holes which can be of different radii £\, • ■ ■ ,e m , define 
Cl £ := n\ U™ x B(ai, Si) and look at 

min / \Du\ 2 (1.7) 

detDu=lJ n£ 

(or minuet Du=i Jq |£>u| n in dimension n), in the limit e — > 0. This also has a reasonable physical 
interpretation: it corresponds to studying the incompressible deformation of a body that contains 
micro-voids which expand under the applied boundary deformation. One may think of the points 
aj as fixed, then they correspond to defects that pre-exist, just as above. Or the model can be seen 
as a fracture model where we postulate that the body will first break around the most energetically 
favorable points aj (see, e.g., the discussion in [1 [SI [691 US HH E3 ED UM \M ED] ) . It can also 
be compared to the core-radius approach in dislocation models [13l \59\ [28] . 

Following the analogy above, we would like to be able to subtract from (]1.7[) a leading order 
term proportional to logj, in order to extract at the next order a "renormalized" term which 
will tell us how cavities "interact" (attract or repel each other), according to their positions and 
shapes. This is more difficult than the problem (11.50 because the condition detDu = 1 is much 
less constraining than |u| = 1. While the maps with |u| = 1 can be parametrized by lifting in the 
form u = e*^, to our knowledge no parametrization of that sort exists for incompressible maps. In 
addition while the only characteristic of a vortex is an integer -its degree-, for incompressible maps, 
the characteristics of a cavity are more complex -they comprise the volume of the cavity and its 
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shape-, and there is no quantization. For these reasons we cannot really hope for something as nice 
and explicit as a complete "renormalized energy" for this toy cavitation model. However we will 
show that we can obtain, in particular in the case of two cavities, some quantitative information 
about the cavities interaction that is reminiscent of the renormalized energy. 

1.2 Method and main results : energy lower bounds 

Our method relies on obtaining general and ansatz-free lower bounds for the energy on the one 
hand, and on the other hand upper bounds via explicit constructions, which match as much as 
possible the lower bounds. This is in the spirit of T-convergence (however we will not prove a 
complete T-convergence result). For simplicity in this section we present the results in dimension 
2, but they carry over in higher dimension. 

To obtain lower bounds we use the "ball construction method", which was introduced in the 
context of Ginzburg-Landau by Jerrard [44] and Sandier [65^ [66] . The crucial estimate for Ginzburg- 
Landau, or more simply S 1 -valued harmonic maps, is the following simple relation, corollary of 
C auchy- S chwar z : 

,9 1 ( f du V d 2 , N 

|£>u| 2 > / — x u = 2vr— (1.8) 

9B(a,r) 2 7Tr \] QB{a,r) ° T J r 

if d is the degree of the map on dB(a, r). Integrating this relation over r ranging from e to 1 yields 
a lower bound for the energy on annuli, with the logarithmic behavior stated above. One sees that 
the equality case in (II. 8ft is achieved when u is exactly radial (which corresponds to u = e in 
polar coordinates), so the least energetically costly vortices are the radial ones. For an arbitrary 
number of vortices the "ball construction" a la Jerrard and Sandier allows to paste together the 
lower bounds obtained on disjoint annuli. Previous constructions for bounded numbers of vortices 
include those of Bethuel-Brezis-Helein [10] and Han-Shafrir [36]. The ball construction method will 
be further described in Section 13.11 

For the cavitation model, there is an analogue to the above calculation, which is also our 
starting point. Assume that u develops a cavity of volume v around a cavitation point a in the 
domain space. By v we really denote the excess of volume created by the cavity (we still refer to 
it as cavity volume), this way the image of the ball B(a,e) contains a volume 7re 2 + v. Using the 
C auchy- S chwar z inequality, we may write 

/ |£>u| 2 > — [ / |Du-r|) . (1.9) 

JdB(a,r) 27rr \JdB{a,r) J 

But then one can observe that Jg B ^ a r \ \Du ■ t\ is exactly the length of the image curve of the circle 
dB(a,r). We may then use the classical isoperimetric inequality 

(Per£(a,r)) 2 > 47r|£(a,r)| (1.10) 

where | • | denotes the volume, and ^(a, r) is the region enclosed by this image curve, which contains 
the cavity, and has volume 7rr 2 + v by incompressibility. Inserting this into (jl.9j) , we are led to 

r |Du| 2 ^ Per 2 (£(a,r)) ^ \E(a,r)\ | ^ 

JdB( a ,r) 2 4vrr r r 
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This is the building block that we will integrate over r and insert into the ball construction, to 
obtain our first lower bound, which is proved in Section [3.11 To state it, we will use the notion of 
weak determinant: 



whose essential features we recall in Section [2.4} as well as Mtiller and Spector's invertibility "con- 
dition INV" [55] which is defined in Section T2.3I (Definition l5l) and which essentially means that the 
deformations of the material, in addition to being one-to-one, cannot create cavities which would 
at the same time be filled by material coming from elsewhere. Even though we have discussed 
dimension 2, we directly state the result in dimension n. 

Proposition 1.1. Let f2 be an open and bounded set in M. n , and £l E = Q\ B(&i,ei) where 
ai,--- , a m E Q, and the £?(aj,£j) are disjoint. Suppose that u £ W 1,n (£l £ , W 1 ) and that condition 
INV is satisfied. Suppose, further, that DetDu = C n in VL £ (where C n is the Lebesgue measure), 
and let Vi := |-E7(aj,£j)| — w n e" (with E(a.i,Ei) as in (jl.lOp ). Then for any R > 



Note that Yli Vi = V is the total cavity volume, which due to incompressibility is completely 
determined by the Dirichlet data, in the case of a displacement boundary value problem. 

Examining the equality cases in the chain of inequalities (jl.9p - (|l.lip already tells us that the 
minimal energy is obtained when "during the ball construction" all circles (at least for r small) are 
mapped into circles and the cavities are spherical. A more careful examination of (j 1.9)) indicates 
that the map should at least locally follow the model cavity map (jl.2p . It is the same argument that 
has been used by Sigalovanathan and Spector \72\ [73] to prove the radial symmetry of minimizers 
for the model with power p < n. 

When there is more than one cavity, and two cavities are close together, we can observe that 
there is a geometric obstruction to all circles "of the ball construction" being mapped into circles. 
This is true for any number of cavities larger than 1; to quantify it is in principle possible but a 
bit inextricable for more than 2. For that reason and for simplicity, we restrict to the case of two 
cavities, and now explain the quantitative point. 

Let ai and a2 be the two cavitation points with |ai — a2| = d, small compared to 1. For simplicity 
of the presentation let us also assume that £i = £2 = £■ The ball construction is very simple in 
such a situation: three disjoint annuli are constructed, B(&i,d/2)\B(a.i,e), B(a2,d/2)\B(&2,e) 
and B (a, R)\B (a, d) , where a is the midpoint of ai and a2 (see Figured]). These annuli can be 
seen as a union of concentric circles centred at ai, a2, a respectively. To achieve the optimality 
condition above, each of these circles would have to be mapped by u into a circle. If this were true, 
the images of B{a\,d/2) and B(a.2,d/2) would be two disjoint balls containing the two cavities, 
call them E\ and E2. By incompressibility, \E\\ = v\ + Tr(d/2) 2 and = v% + ir{d/2) 2 . Then the 
image of -B(a, d) would also have to be a ball, call it E, which contains the disjoint union E\ U E2, 
and by incompressibility 






E\ = vi + V2 + vrd 2 . 



(1.12) 
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Figure 1: Ball construction in the reference configuration 



If d is small compared to V\ and V2 it is easy to check this is geometrically impossible: the radius 
of the ball E\ is certainly bigger than W v\/ir, that of E2 than yjvilit and since E is a ball that 
contains their disjoint union, its radius is at least the sum of the two, hence \E\ > (y/vi + y/v^) 2 - 
This is incompatible with (|1.12p unless 7rd 2 > 1^Jv\Vi- 










II I 













(a) fj, = 0.5 



(b) n = 1.5 



(c) ij, = 3 



Figure 2: Incompressible deformation u : 5(0, d) \ {ai,a2} — > M 2 , d = \&2 — ai|, opening dis- 
torted cavities of volumes v\ + -zre 2 , V2 + 7re|; deformed configuration for increasing values of the 
displacement load (/x : = 



Vl -K<p 2 ) • Choice of parameters: d 



<'_> 
''i 



0.3. 



So in practice, if d is small compared to the volumes, the circles are not all mapped to exact 
circles, the inclusion and disjointness are preserved, but some distortion in the shape of the images 
has to be created either for the "balls before merging" i.e. E\ and E2 - this corresponds to what is 
sketched on Figure [2]- or for the "balls after merging" i.e. E - this corresponds to what is sketched 
in Figure [3] (the situations of Figures [2] and [3] correspond to the test-maps we will use to get energy 
upper bounds, see below Section [L3l . 

A convenient tool to quantify how much these sets differ from balls, which is what we exactly 
mean by "distortion", is the following 
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(a) Reference 
configuration 



(b) Deformed 
configuration, fi = 1 



(c) Deformed 
configuration, fi = 



(d) Deformed 
configuration, = 3 



Figure 3: Incompressible deformation of B(0, d), d := \a. 2 ~ a i|i f° r increasing values of \i :- 
Final cavity volumes v\ and v 2 given by d = 1, ^ = 0.3. 



7rd 2 ' 



Definition 1. T/ie Fraenkel asymmetry of a measurable set E C W 1 is defined as 
D{E) := min J ^ ' — — , with rE such that |-B(x, rjg)| = 

where A denotes the symmetric difference between sets. 

Note that D(E) is a scale-free quantity which depends not on the size of E, but on its shape. 
The following proposition, which we shall prove in Section 13.31 allows to make the observations 
above quantitative in terms of the distortions. 

Proposition 1.2. Let E, E\, and E<i be sets of positive measure in M. n , n > 2 such that E D E1UE2 
and Ei D E2 = , and assume without loss of generality that \E±\ > l-E^I- Then 

TL TL TL 

\E\D(E)~ + |£i|.D(£i)~ + \E 2 \D(E 2 )^= l 



\E\ + \E 1 UE 2 \ 

>C n 



n(n + l) 

\E 2 \ \^ ( (\EA^ + \E 2 \h n -\E\ \ 2( " _1) 



\Ei\ + \E 2 \) \i\Ei\- + \E 2 \^) n - \E X UE 2 \J 
for some constant C n > depending only on n. 

The fact that Ei, E 2 ,E cannot simultaneously be balls is made explicit by the fact that D(E\), 
D(E 2 ), D{E) cannot all vanish unless the right-hand side is negative, which can happen only if \E\ 
is large relative to |-Ei| and l-E^I- The first factor in the estimate degenerates only when one of the 
sets is very small compared to the other. 

Note that such a geometric constraint is also true for more than two merging balls, so in principle 
we could treat (with more effort) the case of more than two cavities, however the estimates would 
degenerate as the number of cavities gets large. 

These estimates on the distortions are useful for us thanks to the following improved isoperi- 
metric inequality, precisely expressed in terms of the Fraenkel asymmetry: 
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Proposition 1.3 (Fusco-Maggi-Pratelli [27]). For every Borel set E C W 1 

A n-l 

PeiE > nu£\E\ — (1 + CD(E)), 

where C is a universal constant. 

In dimension 2, we thus have the improved isoperimetric inequality 

(Per£) 2 > 4ir\E\ +C\E\D 2 (E), (1.13) 
for some universal C > 0. Inserting (I1.13P instead of fjl . 10f> into the basic estimate til . 11 H gives us 

f > |jE(a ' r)l + -\E(a,r)\D 2 (E(a,r)) > - + vrr + -\E(a, r)\D 2 (E(a, r)). (1.14) 

JdB{a,r) ITT V V 

This then allows us to get improved estimates when integrating over r (in a ball construction 
procedure), keeping track of the fact that to achieve equality, all level curves E(a,r) which are 
images of circles during the ball construction would have to be circles. This way, after subtracting 
off the leading order term V{ log ^-, 1 £ we can retrieve a next order "renormalized" term that will 
account for the cavity interaction. This is expressed in the following main result. 

Theorem 1 (Lower bound). Given £1 C M. n a bounded open set, let £l £ := Q \ (B £l (&i) U B e2 (&2)), 
where ai,a2 £ Q, £1,62 > 0, and assume that B £1 (ai) and B £2 (a2) are disjoint and contained in $7. 
Suppose that u S W 1,n (ft £ , W 1 ) satisfies condition INV and Det Du = C n in Q £ . Set 

a = ai + a 2 ^ d:=|ai-a 2 |, vx := \E(&i, ei)\ - u n Z\, v 2 := \E(a 2 , e 2 )\ - u n e%. 

Then, for all R such that B(a, R) C fi, 

R R 
1 ) dx > v x log h v 2 log — 

2ei 2e 2 

- 1 ui n d n \ . I / vi + \ R d 




log min 



V1 + V2 J v\ + v 2 J y \2 n uj n d n J ' d ' maxjex, £2} J 

for some constant C independent ofVt, a\, a 2 , d, v\, v 2 , £\, and £2 (t+ stands for max{0, t}). 

Two main differences appear in this lower bound compared to Proposition ll.il First, the leading 
order term (vi + V2) log ]_ has been improved to v\ log 7- + t>2 log — , which shows that the energy 
goes to infinity as e\ — > or £2 — >■ 0, even if £1 + 62 A 0. This term is optimal since it coincides with 
the leading order term in the upper bound of Theorem [2] below, and in fact it should be possible 
to replace Y2i v i 1°S £ ^2i( Vi F 7 ) m Proposition 11.11 (however, this would require a more 

sophisticated ball construction, and it is not immediately clear how to obtain a general result for 
the case of more than two cavities). Second, and returning to the discusion in dimension two and 
choosing £1 = £2 = £, compared to Proposition 11.11 we have gained the new term 



/ f mm{vi,v 2 }\ 2 nd 2 \ . J Jvi+v 2 R d 



s 



This term is of course worthless unless -^4 — < ( mm i Vl > V2 } \ [ e < mm ^""^ i Under that con- 

V1+V2 \ V1+V2 J — V1+V2 

dition, it expresses an interaction between the two cavities in terms of the distance of the cavitation 
points relative to the data of vi,V2 and e. As v *+ v — > the interaction tends logarithmically to 
+00; this expresses a logarithmic repulsion between the cavities, unless the term log | is the one 
that achieves the min above, which can only happen if logd is comparable to loge. This expresses 
an attraction of the cavities when they are close compared to the puncture scale, which we believe 
means that two cavities thus close would energetically prefer to be merged into one. This suggests 
that three scenarii are energetically possible: 

Scenario (i) the cavities are spherical and the cavitation points are well separated (but not nec- 
essarily the cavities themselves), this is the situation of Figure [3] 

Scenario (ii) the cavitation points are at distance <C 1 but all but one cavity are of very small 
volume and hence "close up" in the limit e — > 

Scenario (iii) "outer circles" (in the ball construction) are mapped into circles and cavities (as 
well as cavitation points) are pushed together to form one equivalent round cavity, this is the 
situation of Figure This seems to correspond to void coalescence (c.f. [8T1 147] ). 

1.3 Method and main results: upper bound 

After obtaining this lower bound, we show that it is close to being optimal (at least in scale). To 
do so we need to construct explicit test maps and evaluate their energy (in terms of the parameters 
of the problem). The main difficulty is that these test maps have to satisfy the incompressibility 
condition outside of the cavitation points, and as we mentioned previously, there is no simple 
parametrization of such incompressible maps. The main known result in that area is the celebrated 
result of Dacorogna and Moser [20J which provides an existence result for incompressible maps 
with compatible boundary conditions. Two methods are proposed in their work, one of them 
constructive, however they are not explicit enough to evaluate the Dirichlet energy of the map. 

The question we address can be phrased in the following way: given a domain with a certain 
number of "round holes" at certain distances from each other, and another domain of same volume, 
with the same number of holes whose volumes are prescribed but whose positions and shapes are 
free; can we find an incompressible map that maps one to the other, and can we estimate its energy 
J |Du| n in terms of the distance of the holes and the cavity volumes? 

We answer positively this question, still in the case of two holes, by using two tools: 

(a) a family of explicitly defined incompressible deformations preserving angles, that we introduce 

(b) the construction of incompressible maps of Riviere and Ye |63} IM] , which is more tractable 
than Dacorogna and Moser to obtain energy estimates. 

We believe it would be of interest to tackle that question in a more general setting: compute 
the minimal Dirichlet energy of an incompressible map between two domains with same volume, 
and the same number of holes, the holes having arbitrary shapes and sizes; and find appropriate 
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geometric parameters to evaluate it as a function of the domains. This question is beyond the scope 
of our paper however and we do not attempt to treat it in that much generality. 



(a) S = 0.1 (b) 8 = 0.4 (c) S = 0.9 

Figure 4: Transition from round to distorted cavities: d = 1, y Vl ^fi. 2 = 1-5, ^ = 0.3. 



Our main result (proved in Section [4. ip is the following. 

Theorem 2. Let ai,a2 G W 1 , v\ > v% > 0, and suppose that d := |ai — a.2\ > £\ + £2- Then, for 
every 5 G [0, 1] i/zere exists a* m £/ie /me segment joining ai and a2, and a piecewise smooth map 
u G C(R™ \ {ai,a2},M") satisfying condition INV, such that DetDu = C n + Vi5 ai + f2^ a2 * n ^ n 
and /or all R > 



Du 



dx < diyi +v 2 + oj n R n ) + i)i ( log — ) + v 2 (log — 



£ 1 / + V £ 2 



(l-5)(log^ +*f»/^log-+ 2 r/^log ° 



B(a«,fi)\(B E1 ( ai )UB E2 (a 2 )) 

4- Hofj), +«o1 I C1 - ^ I lnp- 

d / + V V u i £ i V u i £ 2 

fCi and C2 are universal constants depending only onn). 

If we are not preoccupied with boundary conditions but just wish to build a test configura- 
tion with cavities of prescribed volumes and cavitation points at distance d, then the above result 
suffices. This is obtained by our construction of an explicit family of incompressible maps, which 
contains parameters allowing for all possible cavitation points distances d and cavity volumes V\,V2- 
The feature of this construction is that it allows for our almost optimal estimates, as the shapes 
of the cavities are automatically adjusted to the optimal scenario according to the ratio between 
d, e, y/vi, y/v~2~, their logs, etc, as in the three scenarii of the end of the previous subsection. In 
other words, the construction builds cavities which, when d is comparable to e, are distorted and 
form one equivalent round cavity while the deformation rapidly becomes radially symmetric (as in 
Scenario (hi)); and cavities which are more and more round as d gets large compared to e (as in 
Scenario (i)). For the extreme cases 5 = 1 and 5 = 0, the maps are those that were presented in 
Figures [2] and [3] respectively. The result for intermediate values of 5 is shown in Figure HI 
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/ / 
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\ \ 
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di ai di 




2d 



Figure 5: Geometric construction of domains Oi, 0,2 satisfying 



l«i| 
P2I 



The idea of the construction is the following. Take two intersecting balls B(a.\,p\) and B(a2, P2) 
such that the width of their union is exactly 2d and the width of their intersection is 2d5, and let 
Qi and be as in Figure [5] (the precise definition is given in (|4.2p ). As will be proved in Section 
14.11 for every 5 G [0, 1] there are unique p\ and P2 such that p4 = ^. The cavitation points ai 
and a2 are suitably placed in Q,\ and O2, respectively, in such a way that |ai — a2| = d. It is always 
possible to choose a* between ai and &2 such that J7i U O2 is star-shaped with respect to a*. In 
order to define u in R n \ f^i U Q2 we choose a* as the origin and look for an angle-preserving map 



u(x) = Aa* + /(x) 



X 



x 



A n -1 



vi 



< ! 2 



ISli u 2 | 



By so doing, we can solve the incompressibility equation det Du = 
preserving maps the equation has the same form as in the radial case 



1 explicitly, since for angle- 



det Du(x) : 
which we will see can be solved as 

T(x) 



/ 



n—l. 



X 



9f 
dr 



(x) 



r.n—1 



1. 



X 



X 



.4 



x 



where the function A : S™ -1 — > R is completely determined if we prescribe u on dQ,± U 80,2- Inside 
ill and ^2 the deformation u is defined analogously, taking ai and a2 as the corresponding origins. 
The resulting map creates cavities at ai and a2 with the desired volumes, and with exactly the 
same shape as dQ. 1 and dil.2- For compatibility we impose u(x) = Ax on <90 1 U 80,2- 

In the energy estimate, (1 — 5) log 5 is the excess energy due to the distortion of the 'outer' 
curves u(9B(a*,r)), r G (d, R), and 5 ^ n^f^log + 2 ^/^log j^j is that due to the distortion of 

the curves u(dB(&i,r)), r G (ei,d), i = 1,2 near the cavities. When 6 = 0, fii and VI2 are tangent 
balls, the cavities are spherical, and the second term in the estimate vanishes. The outer curves 
are distorted because their shape depends on that of <9(Qi U CI2), hence a price of the order of 
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(v\ + V2) log 5 is felt in the energy. When 5 = 1, at the opposite end, U O2 is a ball of radius 
d, the deformation is radially symmetric outside f2i U Jl 2 , and no extra price for the outer curves 
is paid. In contrast, the cavities are "D-shaped" (they are copies of <90i and 8^2), and a price of 
order (vi + V2) 2 y/^Tog ~ is obtained as a consequence (in this case the excess energy vanishes as 
^ — > 0, in agreement with the prediction of Theorem [T|). 

Since the last term of the energy estimate is linear in 5, by taking either 5 = or 5 = 1 (and 
assuming R > d) the estimate becomes 

C(vi + u 2 ) mm <^ log — , W — log h 2 A "/ — log — 

[ d \ vi ei V f 1 e 2 

Comparing it against the corresponding term for the lower bound, namel}!!, 

„, . . \ ( v<i \ , R ( v<i \ n ~ x . d 

C(i;i + u 2 ) mm < — log— , — log' 



f 1 / a \v\ ) max{£i,£2/ 

we observe that there are still some qualitative differences. First of all, in the case when e\ <C £2, 
a term of the form log ^ + log ^ is much larger than log max ^ 1 £2 | ■ We believe that the expression 
in the lower bound quantifies more accurately the effect of the distortion of the cavities, and that 
the obstacle for obtaining a comparable expression in the upper bound is that the domains Q± and 
2 m our explicit constructions are required to be star-shaped. For example, in the case d ~ £2, & n 
energy minimizing deformation u would try to create a spherical cavity at ai (so as to prevent a 
term of order log ^- from appearing in the energy due to the distortion of the first cavity), and, at 
the same time, to rapidly become radially symmetric (because of the price of order log 5 due to the 
distortion of the 'outer' circles). Therefore, for values of ite\ <C v\ + V2, the second cavity would be 
of the form B\Bi for some balls B\ and B such that B± C B, \B\\ = v\, and \B\ = V\ + V2- In other 
words, u must create "moon-shaped" cavities, which cannot be obtained if u is angle-preserving. 

In the second place, the interaction term in the lower bound vanishes as ^ — > regardless of 
whether the minimum is achieved at log or at log |, whereas in the upper bound this vanishing 
effect is obtained only for the case of distorted cavities (when log ~ is the smallest). This is 



because when 5 = and v\ 3> V2, the circular sectoio {a* + de l6 ,9 E (§>if")} is mapped to a 
curve Aa* + fiip)^ with polar angles ip ranging almost from to 2tt. This "angular distortion" 
necessarily produces a strict inequality in (jl.9p . so in principle it could be possible to quantify its 
effect in the lower bound. It is not clear, however, whether for a minimizer an interaction term of 
the form (v± + « 2 ) log will always be present (in the case when ^ — > 0), or if the fact that such a 
term appears in the upper bound is a limitation of the method used for the explicit constructions. 

Finally, the factor ^ in front of log ^- and log ^ is raised to a different exponent in each term, 
the reason being that f2i and SI2 play different roles in the upper bound construction. Provided 
5 > 0, when ^ — > the first subdomain is becoming more and more like a circle (its height and its 
width tend to be equal, and the distortion of the first cavity tends to vanish) whereas O2 becomes 



1 When considering boundary conditions, not all values of S can be chosen, see the discussion below. 
2 we assume, e.g., that v± + V2 < 4-kR 2 , in order to illustrate the main point 
3 we state this in two dimensions for simplicity 
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increasingly distorted (the ratio between its height and its width tends to infinity). The factor 
2 \f^ m f ron t °f 1S on ly due to the fact that the effect in the energy of the distortion of the 

cavities also depends on the size of the cavity. 

Dirichlet boundary conditions 

If we want our maps to satisfy specific Dirichlet boundary conditions, then they need to be "com- 
pleted" outside of the ball B(a* , R) of the previous theorem. For that we use the method of Riviere 
and Ye, and show how to obtain explicit Dirichlet energy estimates from it. We consider the 
radially symmetric loading of a ball, but other boundary conditions could also be handled. Let 
a*, 5, pi, p2, Ox, O2 be as before. We are to find R\, R2, and an incompressible diffeomorphism 
u : {Rx < |x - a* I < R 2 } ->■ W 1 such that 

i) Q\ U O2 C B(a*,Ri) and VL\dB(a.*,Ri) coincides with the map of Theorem[2] 

ii) u|aB(a* t R 2 ) is radially symmetric. 































[tffwj 



























(a) S = 0.1, reference configuration, (b) S - 
ir(B% - Rf) = 3.06(vi + v 2 )(l - S) -k{R 2 2 



0.4, reference configuration, 
-Rl)=3.12(vi+v 2 )(l-5) 



(c) S = 0.9, reference configuration, 
■k{RI - R\) = 2.46(ui + v 2 )(l - 5) 






(d) 8 = 0.1, deformed configuration. 
Thick line at u(dB(a*, Ri)). 



(e) S = 0.4, deformed configuration. 
Thick line at u(dB(a*,Ri)) 



(f) 8 = 0.9, deformed configuration. 
Thick line at u(dB(a*, #1)) 



Figure 6: Transition to a radially symmetric map. A larger initial domain is necessary in order to 
create spherical cavities. Parameters: £l = B(0,R2), 



V1+V2 
7T tf 4 



1.5, a =0.3, d 



1, -Ri 



d. 
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Not all values of R\ and R2 are suitable for the existence of a solution, since the reference config- 
uration {R\ < |x — a* I < R2} must contain enough material to fill the space between u(dB(a* , R2)) 
(with shape prescribed by the Dirichlet data) and u(dB(a*, R\)) (whose shape is determined by 
Theorem [21 see Figure [6|). In the case of a radially symmetric loading, the farther Q\ U f2 2 is from 
being a ball, the larger the reference configuration has to be. If 5 = 1 nothing has to be imposed; 
if 5 < 1, we must have that 

u n (R%-R r t)>C(v 1 +v 2 )(l-6) 

for some constant C (see Lemma 14. 5|) . It turns out that the above necessary condition is also 
sufficient, as we show in the following theorem: 

Theorem 3. Suppose that &\, a 2 £ K n and d := |ai — a2| > £1 + £2- Let 5 G [0, 1], v\ > v 2 > 0, 



V s :=2 2n+1 n(v 1 +V2)(l-6), R x > max { [I — , 2d \ , R 2 :={R^ + —. (1.15) 



Un V 0J 



Then there exists a* in the segment joining ai and a 2 and a piecewise smooth homeomorphism 
uG W 1 '°°(W l \ {ai, a 2 }, W 1 ) such that Det Du = C n + v\5 ai + v 2 5 a2 inM n , u\ Rn \ B ( a * R2 j is radially 
symmetric, and for all R > R\ 



1 

n iB(a*,fl)\(B E1 ( ai )UB E2 (a 2 )) 



Du 



R R 

dx < Ci(v\ + v 2 + oj n R n ) + vi log h v 2 log — 

£1 £2 



+ c,^ + ( (1 - „ (>„ g ^ ++t ^± + log i)) . 

The main differences with respect to Theorem [2] are that u is now radially symmetric in 

R n \B(a*,R 2 ) and that log f has been replaced with log a/^JJ = C + log i/^^S^ in the 
interaction term. The proof is presented in Section 14.21 As a consequence we finally obtain 

Corollary 1. Let Q, be a ball of radius R > 2d, with d > E\ + e 2 > 0. Then, for every v± > v 2 > 
i/iere exisi ai, a 2 £ O urai/i |ai — a 2 | = d, and a Lipschitz homeomorphism u : Cl \ {ai,a2} — > W 1 , 
such that Det Du = C n + v\5 ai + f 2 <5a 2 in Q, u\qq = Aid (with X n — 1 := '"^f 2 ), and 



n\(B E1 (a 1 )UB £2 (a 2 )) 



dx < Ci (i>i + u 2 + oo n R n ) + log h t; 2 log — 

£1 £2 



+ C 3 (t* + t») min ((!-*)( log fa + ^X 1 -^) ) + S ( ( & og A + 1 



TOt/i5o :=max|0,l - J- 



The value of <5o is such that <5 > So if and only if u n R n > uj n Ri + V5, with w n i?" := V5 + aj n (2<i) n ; 
the idea is to be able to use Theorem [3] and obtain a final energy estimate depending only on v\ , 
V2, d, ei, £2 and the size |f2| of the domain. 



14 



1.4 Convergence results 

Once we have upper and lower bounds, we are able to show that for "almost-minimizers" one of 
the three scenarii described after Theorem Q] holds in the limit e — > 0. 

Theorem 4. Let Q be an open and bounded set in M. n , n > 2. Let Sj — > be a sequence, that we will 
denote in the sequel simply by e. Let {f2 e } £ be a corresponding sequence of domains of the form Q £ = 
\ UfLi -^e( a i,e)i with m 6 N, ai i£ , . . . , a m]£ £ £1 and e such that the balls £> £ (ai i£ ), . . . , -B £ (a m;£ ) 
are disjoint. Assume that for each i = 1, . . . ,m the sequence {aj i£ } £ is compactly contained in £1. 
Suppose, further, that there exists u £ S W 1,n (£l e , W 1 ) satisfying condition INV, DetZ)u £ = C n in 
Q e , sup e ||u £ || L so(q £ ) < oo and 



Du £ (x) 




dx < ( J>, £ ) log + C ( \n\ + f; vj] , (1.16) 



i=l 



wherdo i>i j£ := |-E(aj j£ , e; u £ )| — w n e n and C is a universal constant. 

Then (extracting a subsequence) the limits = lim £ _>.o a i,e an d v i = hin^o^i,^ i — 1, - - ■ ,771 
are u;eZZ defined, and there exists u G ni< p < n M /1 ' p (0, M. n ) n ^^'"(O \ {ai, . . . , a m },R n ) suc/i i/ioi 

•u £ -um W£? (0 \ {ai , . . . , a m }, R») 

• Det Du £ — Det Du in Q\ {ai, . . . , a m } locally in the sense of measures 

• Det L>u = YlT=i + £ n ^ ^- 
When m = 2, one of the following holds: 

i) «/ai 7^ a 2 and V\,V2 > (assume without l.o.g. v\ > ?; 2y ) ; i/ien 

• t/ie cavities imT(u, ai) and imr(u, a2) (as defined in (12. 3p ) are 6aWs of volume v\, 7J 2 

• |£^(ai )£ , e; u £ ) A imT(u, a,)| — > as e — > /or i = 1, 2 



under the additional assumption that v\ + V2 < 2 n a; n (dist( ai "t a2 , dQ)) n , 



w n |a 2 -ai| n / / o; n (diamQ) 
> Ciexp -C 2 1 H — hlog 




ui + v 2 V V u i + u 2 f l + v 2 

/or some positive constants C\ and C 2 depending only on n; 

ii) i/min{vi,v 2 } = (say t> 2 = 0), then imr(u,ai) (i/ie on/y cavity opened by u) is spherical; 

iii) «/ai = a 2 and t>i,t> 2 > (assume V\ > u 2/ ), i/ien 

• imx(u, ai) is a ball of volume V\ + t> 2 

• |a 2j£ — ai )E | = 0(e) as e — > 



4 Now we write S(aj, e , e; u s ), and not just -E(ai iE ,e), to highlight the dependence on u e . It corresponds to the 
cavity opened by u £ at a^ lE (compare with (|1.10|l and ([273J) ) . 
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• the cavities must be distorted in the following sense (C n being as in Proposition \1.2\) : 

n n n 

liminft , lP ( £ (a,„ e; u g ))-- +i ,, P ( E (a 2 .„s;u.))-- > p ( " (L17) 

In the situation of two cavities, the three cases above correspond to the three scenarii of the 
end of Section [1.21 in the same order. 

The main ingredients for the proof are the comparison of the upper bound (|1.16p with the 
lower bounds Proposition 11.11 and Theorem [H standard compactness arguments, and an argument 
introduced by Struwe [77] in the context of Ginzburg-Landau which allows to deduce from the 
energy bounds sufficient compactness of u e . 

1.5 Additional comments and remarks 

We note first that our analysis works provided that the distance of the cavitation points to the 
boundary does not get small (thus the domain cannot be too thin either). It is an interesting 
question to better understand what happens when they do get close to the boundary, as well as 
the effect of the boundary conditions. 

Second, it follows from our work that it is always necessary to compare quantities in the reference 
configuration with quantities in the deformed configuration, due to the scale-invariance in elasticity. 
For example, we have shown that a large price needs to be paid (in terms of elastic energy) in order 
to open spherical cavities whenever the distance between the cavitation points is small compared 
to the final size of the cavities {uo n d n <C v\ + V2)- If we only know that the cavitation points are 
becoming closer and closer to each other, from this alone we cannot conclude that the cavities 
will interact and that the total elastic energy will go to infinity, as the following argument shows. 
Suppose that u is an incompressible map defined on the unit cube Q C M. n , opening a cavity, and 
satisfying affine boundary conditions of the form u(x) = Ax on dQ, A € R nx?1 . Then, by rescaling 
u and reproducing it periodically, it is possible to construct a sequence of incompressible maps 
creating an increasingly large number of cavities, at cavitation points that are closer and closer to 
each other, in such a way that all the deformations in the sequence have exactly the same elastic 
energy (cf. Ball & Murat [7J; see also [60j H9J [50] ) . This is possible because the cavities themselves 
are also becoming increasingly smaller, with radii decaying at the same rate as the distance between 
neighbouring cavitation points. This example also shows that the strategy of filling the material 
with an arbitrarily large number of small cavities is, in a sense, equivalent to forming a single 
big cavity (there is no interaction between the singularities). Here we complement that result by 
showing that if it is not possible to create an infinite number of cavities, then the interaction effects 
in the energy do become noticeable, and under some circumstances can even be quantified. 

Third, we mention that the idea of partitioning the domain and using angle-preserving maps 
inside the resulting subdomains (as described in Section [L3|) can be used to produce test maps that 
are incompressible and open any prescribed number of cavities (for example by dividing the initial 
domain in angular sectors). 

Finally, we discuss the case p ^ n. It is not clear how to extend the analysis to this case, the 
main reason being that the energy is no longer conformally invariant while the "ball-construction 
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method" is only suited for such cases. To see this in a simple way, let us consider the case of two 
cavities, assuming incompressibility, letting e\ = e 2 - > 0, and let us try to reproduce the steps (|1.8p 
and (II. lip with fll . 14f) . The p-equivalent of (11.14p obtained by Holder's inequality (and by relating 
iDul™- 1 to the area element |(cof Du)u\, see Lemma l3.ip is 



dB(a,r) 



Du(x) 



ve-H*) > Fer(E{a 



(noj n r n x ) 



r \)n-i n-p \e( a. r) » / p 



According to this, when p/nwe may bound from below the energy in B(ai > |)US(a 2) |) (with 



d = |a 2 - ai|) by 

/ 

B(ai,f)UB(a 2 ,f) 



p—n 
1 1 n 



Du(x) 



> [vf- + ^ 



2 



+ C7(«i + ^)S 2 (D(E( ai ,r))^ 
Jo 



1 \ r 



n— p— 1 



where (D(E(a.i,r))) n - 1 ) stands for the average distortion 

(D(E(ai,r)))^) := (vfD(E(a 1 , r))^r + D(£ , (a 2 , r))^ 



Analogously, we can bound the energy in B(a, R) \ B(&, d) (with a = ai + a2 ) by 



/ 



r? 



Dufx) 



> («1 + v 2 ) 



R 



„n-p-l 



+ C(vi + ua) 



£)(£^(a,r))"- 1 r 



i r n-p-l 



B(a,R)\B{a,d) 

and obtain: 



p—n 



B(a,R) 
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Du(x) 



r n-p-l + ( v n + y n _ ^ + 



d^" 



+ C(«i + 



(^(^(a^r)) — )r 



ra— p— 1 



+ / D(E(a,r)) 



— r n-p-l 



III 

Assume that v\ + v 2 is fixed (as is the case in the Dirichlet problem). Let us first consider the 
case p < n. Since the limit e — > is not singular in this case (contrarily to p = n), the problem 
cannot be analyzed by asymptotic analysis. If we guide ourselves only by the second and third 

RE p 

terms (II and III), when p < n we can say the following. The factor v£ + v 2 ~ + v 2 ) n in II 
is minimized when minjui, t^} = 0, hence it motivates the creation of just one cavity (the same 

£ E. 

can be said for the problem with M cavities, because + • • • + vfy is concave and the restriction 
«! + ••• + vm = const, is linear). If the above difference has to be positive, the factor (|)™ 
suggests that the two cavitation points would want to be arbitrarily close, and that the cavities 



17 



will tend to act as a single cavity. This is consistent with the prediction for III; indeed, consider 
the corresponding estimate for p = n: 



-! 


£>u(x) 


n Jfl E nB(a,R) 


y/n- 1 




(D(E(a u r))^) — + f R D(E(a,r))^ — 
r Jd r 

Under a logarithmic cost, it is much more important to minimize the distortions D{E(ai,r)) of the 
circles u(dB(ai, r)), i = 1,2, e < r < ^ near the cavities, rather than the distortion of the outer 
circles D(E(a,r)), r > d. As was discussed before, this leads either to the case of well-separated 
and spherical cavities (scenario (i) in p. [9]), or to the conclusion that if outer circles are mapped 
to circles (scenario (hi)) then the distance between cavitation points must be of order e (Theorem 
I4tiip ). In contrast, When p < n, in the presence of the weight r n ~ p ~ 1 , minimizing the distortions 
D(E(a,r)), r > d gains more relevance compared to the distortion near the cavities. 

For the previous reasons, we believe that the deformations of scenario (i) will not be global 
minimizers, instead the body will prefer to open a single cavity. If multiple cavities have to be 
created, then the cavitation points will try to be close to each other, and the deformation will try 
to rapidly become radially symmetric. The cavities will be distorted and try to act as a single cavity 
(as in scenario (hi), which creates a state of strain potentially leading to fracture by coalescence), 
at distances between the cavitation points that are of order 1 (not of order e). This, in fact, is 
what has been observed numerically [8H Wl\ . 

Let us now turn to p > n. The lower bound reads 



Q e nB(a,R) 



n 



Du(x) 



Vn^l 
+ C(vi + v 2 )« 



r ( V , -I- Vo ) „ V P f 2 , p 

dx + 1 1 21 R n ~ p > (v? + ) / r n -v~ x + {y x + v 2 ) » d n ~P 

Je ^~ 



p — n 



li 



(D{E(a i ,r))^)r n - p - 1 + / D(E(a,r))^r n - p - 1 

d 



R 



This time the limit e — > is singular, even more so than for p = n. The factor v" + is now 
minimized when the cavities have equal volumes. Regarding d, the hrst term prefers small distances 
(d = 2s) while the second prefers d — > oo; since (v\ + v 2 ) n > v" + , it can be said that II has 
a stronger influence, hence d large should be preferred^! With respect to the third term, it is now 
much more vital to create spherical cavities (so as to minimize the first of the two integrals) than 
when p = n. This implies that it is scenario (i), rather than (ii) or (iii), which should be observed. 

The case p < n, therefore, should favour a single cavity and coalescence, p > n should favour 
many cavities and splitting, and both situations are possible in the borderline case that we have 
studied: p = n. 



5 although in order to be sure it would be necessary to compute the energy in the transition region B(a,d) \ 
(B(a 1 ,f)uB(a 2 ,|)) 
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1.6 Plan of the paper 

In Section [2] we describe our notation and recall the notions of perimeter, reduced boundary, 
topological image, distributional determinant, and the invertibility condition INV. In Section [3] 
we begin by extending (|1.14p to the case of an arbitrary power p and space dimension n (Lemma 
13. ip . In Section [3.11 we prove the lower bound for an arbitrary number of cavities using the ball 
construction method (Proposition 11.1) ). In Section [3.21 we prove the main lower bound (Theorem 
[I]) and postpone the proof of our estimate on the distortions (Proposition II .2|) to Section 13.31 The 
energy estimates for the angle-preserving ansatz are presented in Section \A. II and proved in Section 
14.31 In Section H2] we show how to complete the maps away from the cavitation points so as to fulfil 
the boundary conditions, and in Section 14.41 we comment briefly on the numerical computations 
presented in this paper based on the constructive method of Dacorogna &; Moser [20J. Finally, the 
proof of the main compactness result and of the fact that in the limit only one of the three scenarii 
holds (Theorem HJ) is given in Section [5j 
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2 Notation and preliminaries 

2.1 General notation 

Let n denote the space dimension. Vector-valued and matrix-valued quantities will be written in 
bold face. The set of unit vectors in R n is denoted by S n-1 . Given a set E C W l , A > and 
h € M. n , we define XE := {Ax : x G E} and + h := {x + h : xE E}. The interior and the 
closure of E are denoted by Int E and E, and the symmetric difference of two sets E± and E2 by 
E1AE2. If Ei is compactly contained in E2, we write E\ CC Ei- The notations B(x, R), B^(x) are 
used for the open ball of radius R centred at x, and B(a,R), .Br (a) for the corresponding closed 
ball. The distance from a point x to a set E is denoted by dist(x, E), the distance between sets by 
dist(-E?i, E2), and the diameter of a set by di&mE. 

Given A an n x n matrix, A T will be its transpose, det A its determinant, and cof A its 
cofactor matrix (defined by A T cof A = (det A)l, where 1 stands for the n x n identity matrix). 
The adjugate matrix of A is adj A = (cof A) T . 

The Lebesgue and the /e-dimensional Hausdorff measure are denoted by C n and 'H k , respectively. 
If E is a measurable set, £ n (E) is also written \E\ (as well as |/| for the length of an interval I). 
The measure of the £>dimensional unit ball is uik (accordingly, % n_1 (<9B(x, r)) = n^r™" 1 ). The 
exterior product of 1 < k < n vectors ai, . . . , a& S M n is denoted by ai A • • • A a^ or f\ i=1 a^. It is 
/c-linear, antisymmetric, and such that |ai A • • • A a&| is the fc-dimensional measure of the fc-prism 
formed by ai, . . . a fc (see, e.g., [Ml ESI [331 E])- In particular, |x| 2 = |x • e| 2 + |x A e| 2 for all xel" 
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and e 6 S n_1 . With a slight abuse of notation, when k = n the expression ai A ■ ■ ■ A a n is used to 
denote the determinant (in the standard basis) of the matrix with column vectors ai, . . . , a n G R n . 

The characteristic function of a set E is referred to as xe, and the restriction of u to E as u\e- 
The sign function sgn : M —> { — 1, 0, 1} is given by sgnx = x/\x\ if x ^ 0, sgnO = 0. The notation 
id is used for the identity function id(x) = x. The symbol j E f stands for the integral average 
j]kf E f- The support of a function / is represented by spt /. 

The space of infinitely differentiable functions with compact support is denoted by C£°(£l), and 
the LP norm of a function / by ||/||lp. Sobolev spaces are denoted by W 1 'P(fi,R n ), as usual. The 
Hilbert space W 1,2 (Q, M 71 ) is denoted by i7 1 (0, IR n ). The weak derivative (the linear transformation) 
of a map u 6 W 1 ' p (£l,W l ) at a point x G W 1 is identified with the gradient Z)u(x) (the matrix of 
weak partial derivatives). 

Use will be made of the coarea formula (see, e.g., \24\ [23"1 Q]): if E C M. n is measurable and 
<j> : E ->-M is Lipschitz, then for all / G L l {E) 

[ /(x)|Zty(x)|dx= I™ if /(xJd^W] di. 

je j-oo yj{xeE:0(x)=<} y 

2.2 Perimeter and reduced boundary 

Definition 2. The perimeter of a measurable set E C W 1 is defined as 
Per E := sup j ^ div g(y) dy : g G C\ (R' 



mn ),llKlU < 1 



Definition 3. Given yo G M n and a non-zero vector v G M n , we define 

H + (y , v) := {y G R n : (y - y ) • > 0}, iT (y 0) i/) := {y G M n : (y — yo) • v < 0}. 

T/ie reduced boundary of a measurable set E C W 1 , denoted by d*E, is defined as the set of points 
y G W 1 for which there exists a unit vector v G W 1 such that 

l.„ |gnH-(y y )nfl(y.r)| = l ^ lml 1^ n H*(y „) n B(y,r)\ = ft 
r^o+ l-B(y,OI 2 |S(y, r)| 

If y G i/ien 1/ is uniquely determined and is called the unit outward normal to E. 

The definition of perimeter coincides precisely with the % n-1 -measure of the reduced boundary, 
as follows from the well-known results of Federer, Fleming and De Giorgi (see, e.g., |24 l l83 t [23"| , [T]l 6 l. 

2.3 Degree and topological image 

We begin by recalling the notion of topological degree for maps u that are only weakly differentiable 
[561 [2S1 HI QH- 

If u G TU 1 'P(Q,lR ri ) and x G W 1 , then, for a.e. r G (0, oo) with 8B(x,r) C 



6 When Per E = oo, the result is true if we consider the measure-theoretic boundary, as defined in [231 Th. 5.11.1]. 
For sets of finite perimeter the two notions of boundary coincide H n ~ -a.e., thanks to a result of Federer [24] (also 
available in [U Th. 3.61], [23l Lemma 5.8.1], or [HI Sect. 5.6]). 
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(Rl) u(z) and Du(z) are denned at W 1 1 -a.e. z G dB(x,r) 
(R2) u| efl(x>p) G W 1 -P(95(x ) r) > R") 

(R3) D(u|g B ( xr ))(z) = (-Du(z))| Tz (g B ( xr )) (the n-dimensional and the tangential weak derivatives 
coincide; T z (dB(x,r)) denotes the tangent plane) for "H n_1 -a.e. z G dB(x,r) 

(this follows by approximating by C°° maps and using the coarea formula). If, moreover, p > n— 1, 
then, by Morrey's inequality, there exists a unique map u G C°(dB(x.,r)) that coincides with 
u las(x,r) H n -&.e. With an abuse of notation we write u(9B(x,r)) to denote u(dB(x, r)). 

If p > ra — 1 and (R2) is satisfied, for every y G ]R n \ u(dB(x, r)) we define deg(u, dB(x, r),y) as 
the classical Brouwer degree [68j ES] of u laB(x,r) with respect to y. The degree deg(u, dB(x, r), ■) 
is the only L 1 (W l ) map [56l E2] such that 

[ deg(u,aS(x,r),y)divg(y)dy = f g(u(z)) ■ (cof Du(z))i/(as) dW^^z) (2.1) 

JB™ JdB{x,r) 

for every g £ C 1 (lR n ,lR n ), f(z) being the outward unit normal to dB(x,r). 

For a map u G W 1,p (£l, M. n ) that is invertible, orientation-preserving, and regular except for 
the creation of a finite number of cavities, deg(u, dB(x, r), y) is equal to 1, roughly speaking, only 
at those points y enclosed by u(dB(x,r)). Because of this, the degree is useful for the study of 
cavitation, since we can detect a cavity by looking at the set of points where the degree is 1, but 
which do not belong to the image of u (they are not part of the deformed body). This gave rise to 
Sverak's notion of topological image |79j . 

Definition 4. Let u G W 1 ' p (dB(x,r),W l ) for some x G K ? \ r > 0, and p > n - 1. Then 

im T (u,fl(x,r)) :={y el" : deg(u, dB(x, r), y) 0}. 

It was pointed out by Miiller-Spector [55J Sect. 11] that Sobolev maps may create cavities in 
some part of the body, and subsequently fill them with material from somewhere else (even if they 
are one-to-one a.e. [3]). In order to avoid this pathological behaviour, they defined a stronger 
invertibility condition, based on the topological imaged. 

Definition 5. Let u G W 1,P (Q,, W 1 ) with p > n — 1. We say that u satisfies condition INV if 

i) u(z) G imr(u, -B(x, r)) for a.e. z G B(x, r) 

ii) u(z) G W 1 \ im T (u, B(x, r)) for a.e. z G \ B(x,r) 

for every x£l™ and a.e. r G (0, oo) such that u|g B / Xjr ) G W l,p (dB(x,r),W n ). 

7 The original definition of condition INV in [551 Sect. 3] required that i) and ii) were satisfied only for a.e. r 6 (0, oo) 
such that _B(x,r) C f2. Here we impose i) and ii) for a.e. r £ (0, oo) such that 9_B(x,r) C As explained in [37], 
this modification is necessary when considering perforated domains, due to Sivaloganathan & Spector's example of 
leakage between cavities |74l Sect. 6]. 
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In the following proposition we summarize some of the main virtues of condition INV. We add 
a sketch of the proof to make it easier for the interested reader to compile the different ideas and 
conciliate the different notation in [79], [551 Lemmas 2.5, 3.5 and 7.3], [17^ Lemmas 3.8 and 3.10], 
[391 Lemma 2], and \40\ Prop. 6 and Lemma 15]. 

Proposition 2.1. Let u G W 1,P (Q, M n ) with p > n — 1 satisfy det Du > a.e. and condition INV. 
Then, for every x £ W 1 there exists a full- C 1 -measure subset R x of {r G (0, oo) : dB(x, r) C £1} for 
which (R1)-(R3), conditions^ \v§) of Definition^ and the following properties are satisfied: 

i) deg(u, dB(x, r), y) G {0, 1} for every y G W 1 \ u(dB(x, r)) 
ii) d* imT(u, -B(x, r)) = u(dB(x,r)) up to "H n_1 -null sets 

Hi) Per (im T (u,£(x,r))) = / |(cof Du(z))i/(z)|d^ n_1 (z) 

JdB{*,r) 

iv) |im T (u,B(x,r))| = - f u(z) • (cof Du(z))i/(z) dft n_1 (z). 

Moreover, for every x, x' G M n and every r G i? x , r ' £ R*' 

v) im T (u,B(x,r)) C im T (u, B(x', r')) if B(x,r) C B(x',r') 

vi) im T (u, B(x, r)) n im T (u, 5(x', r')) = if B(x, r) n B(x',r') = 0. 

Proof. Call S7q the set of x G SI for which there exist w G C 1 (IR n ,M n ) and a compact set K C £1 
such that 

\K n -Bfx, r)| 

lim — — — = 1, u\ K = w\k, and Du\ K = Dw\ K . (2.2) 

r->0+ |-d(x, r)| 

Since u G VF 1,p (f2, M n ), it is possible to find (combining Federer's approximation of approximately 
differentiable maps by Lipschitz functions, Rademacher's theorem, and Whitney's extension the- 
orem, see, e.g., [231 Cor. 6.6.3.2], [SJ Thms. 3.1.8 and 3.1.16], [55j Prop. 2.4], [39J Lemma 
1]) an increasing sequence of compact sets {Kj}j £ ^ contained in fi, and a sequence {wj}j g ^ 
of maps in C 1 (IR n ,lR n ), such that = w j|ifji Vuj|^ = Dw\Kj, and \Q \ Kj\ < X for 

each j G N. By Lebesgue's differentiation theorem, \Kj \ K'-\ = where K'- := {x G Kj : 
lim r .^ +(r- n | J B(x, r) \ K\) = 0}. Since ft => UjeN^j' it follows that \Q \ O | = 0. 

Define i? x as the subset of {r G (0, oo) : dB(x,r) C fi} for which (R1)-(R3), conditions |T])-|TiJ) of 
Definition [H and the following properties are satisfied: 

(R4) ^"- 1 (aS(x,r)\O ) =0 

(R5) det£>u(z) > for % n-1 -a.e. z G dB(x,r). 

The fact that \{r G (0, oo) : dB(x, r) C £1} \ R x \ = is a consequence of the coarea formula and of 
the discussion before Definition [U For this choice of i? x we have that the properties listed in the 
proposition are satisfied for all (not only for a.e.) r G i? x . This follows from (|2.ip . the fact that 
U|n is one to one (by [551 Lemmas 3.4 and 2.5]; only minor modifications are required, see [391 
Lemma 2] if necessary), and a careful inspection of the proofs of |55l Lemmas 2.5, 3.5 and 7.3]. □ 
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By Proposition 12.11^]) the topological image of B(x,r) can be defined for all x E M. n and all 
r > such that {z : r < \z\ < r + 5} C f2 for some 5 > (not only for radii r E -R x )- Indeed, since 
the sequence {imx(u, S(x, r)) : r E i? x } is increasing for every x E M n , we may define 

£(x,r):= P| im T (u,£(x,r)). (2.3) 

r'>r 
r'Gflx 

Whenever explicit mention of u is necessary (such as in Theorem 2] where sequences of deformations 
are considered), we write E(a, r;u). Finally, if a point a E M. n is such that D(a, 5) \ {a} C O for 
some 5 > 0, we define its topological image as E(ai,0), and denote it by imx(u, a). 

2.4 The distributional determinant 

It is well known that the Jacobian determinant of a C 2 vector- valued map u : Q C W 1 — > W 1 has a 
divergence structure. When n = 2 or n = 3, this is 

det £>u = ui ; iu 2 ,2 - u 2 ,iui i2 = ("i"2,2),i - ("i"2,l),2 



det Du = itx i 



"l 



"2,2 "2,3 

""3,2 "3,3 

"2,2 "2,3 

"3,2 "3,3 



+ U1.2 



"2,3 "2,1 
"3,3 "3,1 



+ "1,3 



"2,1 "2,2 
"3,1 "3,2 



+ "1 



"2,3 "2,1 
"3,3 "3,1 



+ "1 



"2,1 "2,2 
"3,1 "3,2 



.3 



where Uij denotes the j-th partial derivative of the i-th component of u. In higher dimensions, we 
may write det-Du = Div((adj Du)^). 

One of the main ideas in Ball's theory for nonlinear elasticity [2] is that if the divergence is taken 
in the sense of distributions, the right-hand side of the above expressions is well defined for maps 
that are only weakly differentiable. This motivated his definition of the distributional determinant 
of a map u E W 1,n-1 (n, M n ) n R n ) as the distribution Det Du E D'(fi) given by 



(DetDu,0) := 



1 

" J si 



u(x) • (cof Du(x))D(f>(x) dx, 



E C c °°(tt) 



(2.4) 



(see also [HU [161 QUI EH [2IJ [TT] and references therein for subsequent developments and for the role 
of Det-Du in compensated compactness, homogenization, liquid crystals, and superconductivity). 
If a map u E W ,p (f2,R n ), p > n — 1, satisfies condition INV, then u(z) is contained in the 
region enclosed by u(dB(x, r)) for every x E M n , a.e. z £ On B(x,r), and a.e. r > such that 
dB(x,r) C 0. Consequently, u E LJ^ (J), R n ), and the distributional determinant is well defined. 



Proposition 2.2 (cf. [55], Lemma 8.1). Let u E W lj, (n,R n ), p > n 
and condition INV. Then 

i) Det Du = (det Du)C n + where fj, s is singular with respect to C n 

ii) \E(x, r) \ imx(u, -B(x, r))\ = for every x E K n anrf r E i? x 



1, satisfy det-Du > a.e. 
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Hi) \E(x, r 2 ) \E(x., ri)| = Det Du(A rit1 . 2 ) for all r\ > and r 2 > such that the annulus A rir2 := 
{x G ]R n : ri < |x| < r 2 } is contained in Q. 

Proo/. Let x£l" and set 5:={r£ (0, oo) : <9P(x, r) C 0}. The map 

cj( r ) :=- [ u(z) • (cof Ou(z))i/(z) dW n_1 (z), r G P x 

n JdB{x,r) 

belongs to L l (S). Suppose [ri,r2] C S for some r\,r 2 G P x . For 5 > let <^(z) := V<5(l z ~~ x I)j 
where ips G C^°([0, oo)) is such that ipg = 1 in (ri + 5, r2 — (5), ^ = in [0, r\\ U [7-2,00), and 
^ll^lloo < 2. It is clear that 4>s — > XA n , r2 pointwise as 5 — > + , and that 

(Det Du, fa) = co(r 2 ) - w(n) + ^(r)(w(n) - w(r)) + <ty£(r)(w(r 2 ) - w(r)). 

The proof follows from |55| Lemma 8.1], Proposition 12. Hiv|) -fvj) . and Lebesgue's differentiation the- 
orem aplied tow. □ 



3 Lower bounds 

The following is the basic estimate that allows us to relate the elastic energy to the volume and 
distortion of the cavities. It extends (|1 . 14[) to an arbitrary exponent p and dimension n. 

Lemma 3.1. Suppose that u £ W 1,P (Q, W 1 ), p > n — 1, satisfies detDu > a. e. and condition 
INV. Then, for every x G £1 and r G Pt x ( defined in Proposition \2.1\) , 



i 



9B(x,r) 



f(u|8fl(x,r))(x) 



dU 



n-1 



X > 



\E(x,r) 



[l + CD(E(B(jc,r))) n - 1 . 



|P(x,r) 

Equality is attained only if u\dB(x,r) * s radially symmetric. 

Proof. Given x G K n , r > and z G <9P(x, r) such that Du(z) is well defined, we have that 
|(cof£>u(z))i/(z)| = |(Du(z))ei A • • • A (Du(z))e„_i| < |(Du)ei| • • • |(P>u)e n _i| 

2 r 



< (n-l)^r (|(P>u)ei| 2 H h \ {Du)e Ti 



■"n—\ 



{ei, . . . , e n _i, f(z)} being an orthonormal basis of R n with v(z) := (z — x)/r. Equality holds only 
if |(Du)ej| = |(P>u)ej| and (Du)ej _L (Du)ej for i / j, as in Sivaloganathan-Spector [721 [73]. If 
r G P x , by Propositions 12. HtiT|) . I2.2iij) . and 11.31 we obtain 



£ 



8B(x,r) 



D(u 



dB(x,r) 



n-1 



The conclusion follows by Jensen's inequality. 



□ 



24 



3.1 Ball constructions, the case of multiple cavities 

In this Section we prove Proposition 11.11 (our first lower bound, valid for an arbitrary number of 
cavities). We start by introducing the necessary notation, and by recalling the ball construction 
method in Ginzburg-Landau theory, following the presentation in |67j . 

Collections of balls will be denoted by expressions with B. If B is a ball, r{B) denotes its radius. 
If B is a collection of balls, then r(B) = £ BeB r(£). If A > 0, XB := {XB : B G £>}. We use \JB 
to denote the union UbsB B oi a collection of balls. Given a measurable set A and a collection of 
balls B, we denote {B n A : B G B} by A n B. Given F : W 1 x (0, oo) R, we regard J 7 as a 
function defined on the set of all balls (cf. [671 Def. 4.1]), and write J~{B) for J-(x., r) if B = B(x, r) 
(or £?(x, r)). Also, we write J 7 ^) for XlseB -^( B ) if S is a collection of balls. 

Proposition 3.2 (cf. [67], Th. 4.2). Zet So 6e a finite collection of disjoint closed balls and let 
to := r(Bo). There exists a family {B(t) : t > to} of collections of disjoint closed balls such that 
B(to) = Bq and 

i) For every s>t>t ,{J B{t) C |J B(s). 

ii) There exists a finite set T such that if\t\,t<2\ C [to,oo) \ T, then Bfo) = j^B(ti). 
Hi) r{B{t)) = t for every t > to. 

We point out that we chose a different parametrization from the one in \67\ Th. 4.2]. Here t 
corresponds to e t there. 

Definition 6 ([67], Def. 4.1). We say that a function T : W 1 x (0, oo) — > M is monotonic (when 
regarded as a function defined in the set of balls) if J-(x,r) is continuous with respect to r and 
J-(B) < J-(B') for any families of disjoint closed balls B,B' such that [JB C \JB'. 

Proposition 3.3 (cf. [67], Prop. 4.1). Let T : W 1 x (0,oo) — > M be monotonic in the sense of 
Definition® Let Bq and {B(t) : t > to} satisfy the conditions of Proposition HOI Then, 

fS px-p Af 

HB{s))-F{B )> Y, r ^^~T C 3 - 1 ) 



<0 B(x,r)eB(t) 



for every s > to, and for every B G B(s) 



fS AT Af 

T{B)-T{B Q C\B)> Y r ^( x ' r )y- ^ 

Jt ° B(x,r)eB{t)nB r 

Lemma 13.11 applied to "^"(x, r) = f B ^ x ^ ^ ^==[ | — 1^ dx and Proposition 13.31 immediately 
imply the following result (stated without proof). 

Proposition 3.4. Suppose that u £ W 1 ' p (ft,W n ) with p > n — 1 satisfies detDu > a.e. and 
condition INV. Suppose, further, that Bq and {B(t) : t > to} satisfy the conditions of Proposition 
3.2\ Then, for every s > to such that £l s := |J B(s) \ |J Bo C fi, 

Dufx) 



1 

n Jn 3 



I E \ B \ l^rO- + CD{E B ))£i 
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where Eb denotes E(x,r) for B = B(x,r). Analogously, for every B £ B(s) 



1 

n 



BMJBi 



Du(x) p \ 

1 dx > 



E 1^1 

*° B>eB(t)nB \\ B \ n ) 1 



Proposition 11.11 finally follows from Proposition 13.41 and the incompressibility constraint: 

Proof of Proposition \l.l[ Let A := {i : B(a.i,R) C 0}, to := r (Bo) = SieA 6 *' anc l := 
\J ieA B ei (a.i). Let {S(t) : t > to} be the family obtained by applying Proposition 13.21 to Bo. 
Then, applying Proposition 13.41 if U^( s ) c ^> 



- f ( ^4 " - l) dx > f £ ((|£ B | - + C|£b|0(£b) 



dt 



(3.3) 



BeB(t) 

Proceeding as in the proof of Proposition 12.21 and using incompressibility we obtain 



E B \ U E{ai,£i 

SLitB 



DetDu [ B\ \J B £i (ai) J = \B\ - E cu n e™, 



hence, by the definition of V{ in the statement of the proposition, 



\E B \ - \B\ 



a; SB 



E = E 



(3.4) 



Combining ()3.3[) and (|3.4p we obtain 



UBW 



Du(x) 



y/n=T 



E ^ 

>i,S( ai ,R)cn E 



logf + C [ S I E |^s|I>(£7fl)^ ] 
t0 Jt ° \Bes(t) / 



dt 

T' 



Let so := sup{s £ [io>-R) : U^(' s ) C O}. If so = -R, the claim is proved. Otherwise, from 
Proposition 13.21 we deduce that there exists a ball B(a,r) £ i3(so), of radius r < so, containing at 
least one aj, i £ A, such that B(a,r) n <9f2 ^ 0. The proof is completed by observing that 



R < dist(aj, 90) < |aj — a| +dist(a, c?fi) < 2sq. 



□ 



3.2 The case of two cavities: proof of Theorem [T] 

In this section, we prove Theorem [T] assuming Proposition II. 2\ whose proof is postponed to Section 



We will need the following lemma. 
Lemma 3.5 (Modulus of continuity of the distortion). Let E,E' C M n be measurable. Then 
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i) \\E\D{E) - \E'\D{E')\ < 2\EAE'\ 

ii) \E\D{E)— -\E'\D(E')— < 2^^\EAE'\. 

Proof. Let B' be a ball such that \B'\ = \E'\ and \E'\D(E') = \E'AB'\. For all measurable sets B 

\EAB\ - \E'\D(E') = \\xe - Xsllii - \\xe> - Xfi'11-L 1 < llxs ~ XK'llz 1 + Wxb ~ Xb>\\ia- 
Testing with concentric balls, and taking the minimum over all balls B with |5| = yields 

\E\D(E) - \E'\D{E') < \\xe ~ XE'WlA + ll^l " l^'H 

(||X-B — Xb'IIl 1 = 11-^1 ~~ l-^'ll since B and B' are concentric). Combining this with the fact that 
\\E\ - \E'\\ = HIxeIIli - I|Xb'||li| < IIXE-XE'llii, weobtainE]). 

Property |n|) follows froinG}, the mean value theorem, and the fact that D(E) < 2 for all E (a 
direct consequence of its definition). To be more precise, suppose that \E\ > \E'\, then 

\E\D(E)^ - \E'\D(E')^ 

1 n 

\E\~—{\E\D(E))— - \E'\- — (\E'\D(E')) — 



< \E\- 
2n 

< 



(\E\D(E))& - {\E'\D(E'))^ +{\E'\D{E'))^ \E\ 



\E'\ 



n — 1 

completing the proof. 



\E\- — (max{\E\D(E), \E'\D(E')})^\EAE'\ + 



2~ 

n — 1 



\\E\-\E% 



□ 



We now proceed to the proof of Theorem [TJ As in (13. 4h , by Proposition 12.21 we have that 
\E(B)\ = \B\ + Y^i- ai eB v i f° r au balls B with 3i3 C f2 e . Hence, Lemma 13 . 1 1 implies that 



n 



dB(x,r) 



Du(x) 



1 dH n - L (m) > 



E 



. i: aiS_B(x,r) 



n \ 1 

+ C , |£(x,r)| J D(£(x,r))- 1 ] - (3.5) 



for all x£R" and all r e R x . Given R > d such that B(a, R) C fi, let 



:= % 2 (ai) \ B ei (ai), A 2 := 5 d/2 (a 2 ) \ B £2 (a 2 ), A 3 := Sfl(a) \ B d (a). 

By considering that £l e n -B(a, i?) D A\ U A 2 U A3 and integrating successively in each annulus, we 
obtain 



1 



n 



n e nB(a,R) 



Du(x) 



d d . R 

1 dx > vx log h v 2 log h (v 1 + u 2 ) log — 

2ei 2e 2 a 



(3.6) 



r rf / 2 « dr Z^/ 2 « dr 

+ C \E( ai ,r)\D(E( ai ,r))^— + C \E(a. 2 ,r)\D(E(a 2 ,r))*=i — 

Jei T J 62 T 



n A f 

+ C \E( a ,r)\D(E( a ,r)) — — . 

J d ' 
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Proposition 11.21 applied to E\ = E(&i, |), E 2 = E(a 2 , f), and E = E(a,r), r G (d, R) gives 

n 

\E(a,r)\D(E(a,r)) n - 1 



n(n+l) 
2(n-l) 



>_ c( „ + „ 2) ^L±^M — ( Wl;ffl> ) - (3 . 7) 

n n 

- \E(a u d/2)|£>(£7(ai,d/2)) - |£(a 2 , d/2)|£>(£(a 2 , d/2)) . 

i i 



Define 5(/?i,/? 2 ) := (Pi + T - (Pi + AO (when n = 2, 5 (/3i,/3 2 ) = 2^A^). Using that 



Ei\ = v i + i i = 1,2 we may write 



ii uj d n 

(\Ei\n + |£ 2 |nf = 5(1^1, |E 2 |) + (|£?x| + \E 2 \) = 5(1^11,1^21) + 2 • ~^ r +v 1 +v 2 . (3.8) 



Estimate (|3.7p is meaning ful if \E(a,r)\ < (|#i|™ + \E 2 \n) n , i.e. if 



uj n d n (jj n d n \ 0J n d n 
^ < u n r n < g ( V! + -^,v 2 + J + ^zr (3-9) 



(since 5 is increasing in /?i and /? 2 and g(/3, /?) = (2 n — 2)/3, the inequality holds at least for r = d). 
Define p as the radius for which uj n r n is in the middle of the two extremes in (|3.9D . 



:= (2«- + 1)^ + ±^ 1 + ^, U2 + ^| (3.10) 

For all r £ (d, mm{p, R}) we have that E(a,r) C E(a,p), hence 

1 1 ) rl n 

\E(a, r)\ < u n p n + v 1 + v 2 = -g(\Ei\, \E 2 \) + (2™" 1 + 1)^— + Vl + y 2 . (3.11) 

Noticing that g is 1-homogeneous, combining (|3.8p and (13. llf) we obtain 

(l^l^ + l^l^r-l^r)! > i 5 (|^i|,|^ 2 |)-(2"- 1 + l-2)^ _1 2"" 1 -! 



d^+i^)»-i^u^r 5(1^1,1^1) 2 

' ,J , , 1 ,J by 



n-1 _i n— 1 y \ 

5(1 + x, 1 + y) = £ Q ((1 + x) k (l + y) n - k ) ~ > £ (JJ (1 + fca:)£(l + (n - %)* 

n— 1 / s 

^ E U (1 + x) " (1 + y) " - {r " 2)(1 + x + y) " (3 - 12) 
( with x = - 1 = & and y = &) to obtain 



(|£?i|n +|E 2 |n)"-|^(a,r)| 
{\E^ + \E 2 \l) n - \E X VJE 2 \ 



n(n + l) 

n(n + l) / \ 2(n-l) 

1 2 n_1 - 1 1 "("+!) 



2(n-l) 

> 



V 2 (2»-2)(l + 2»ad£)"_ 



> 4 ^"-i) 
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On the other hand, \E\ U E 2 \ < 2(vi + v 2 ) (because oj n d n < v\ + v 2 ), and since \E±\ > v\ and 
\E 2 \ — v 2, we can substitute with m1 ^^' 2 ^ in (|3.7|) . Hence, for all r G (cZ, rrrinjp, i?}), 

all si G (ei,d/2) and all S2 G (e 2 ,d/2), 

\E( a ,r)\D(E(a,r))^ + |£(ai, ai)|D(^( ai> si))^ + |£(a 2 , s 2 )|£>(£(ai, ai)) ^ (3.13) 
> C(n)( Vl + U2 ) ( m ^^ 2} ) " E |l^( a - s i)P(^(ai, s t ))^ - \E(ai, $)\D{E(ai, |))^ 



Denoting E^a^si), £'(a 2 ,s 2 ), and E(a,r) by i£ Sl , E? S2 , and E r , from (|3,6p we obtain 



1 

n Jn £ nB(a,R) 



Du(x) 



1 ) dx > ui log h i>2 log — (3.14) 

1E\ l£ 2 



VrV^T 

+ C inf (|£ r |£(£ r )^ + 

re(d,min{p,R}) V / {d d £ ) 

Sie(£i,d/2) 

with e = max{ei,e 2 }. In order to estimate log ^, from ()3.10p and (|3.12p we find that 

P -> O— (n+l)„ / 1 1 2 n t)i 1 1 2 n t' 2 \ \ (n— 1 _ n-n\ / -i , nnt>i+2^.\ 71 •> /I — 9 1 ~ nN l f ^1+^2 A n 

The proof is completed by combining (|3.13j) and (|3.14p with Lemma 13.51 

3.3 Estimate on the distortions 

This section is devoted to the proof of Proposition 11.21 

Lemma 3.6. Let q > 1 and suppose that E, E\, and E 2 are sets of positive measure such that 
E D Ei U E 2 and E X C\E 2 = 0. Then 

\E\D{Ef + mDjE^ + \E 2 \D(E 2 )i > / || XB - XBl - XJ?2 || Ll - (LBj - \B X \ - \B 2 \) \ q 

\E\ + |£iU£ 2 | ~~ b,B\,b% \ \E\ + \E X UE 2 \ J ' 

where the minimum is taken over all balls B, B\, B 2 with \B\ = \E\, \B\\ = \E\\, \B 2 \ = \E 2 \. 

Proof. Let B, B±, B 2 attain the minimum in the definition of D(E), D(Ei), D(E 2 ), that is, suppose 
that |J3| = \E\, \Bi\ = \Ex\, \B 2 \ = \E 2 \ and 

\E\D(E) = \EAB\, \E l \D(E 1 ) = \E 1 AB 1 \, \E 2 \D(E 2 ) = \E 2 AB 2 \. 

Since \b - Xb, - Xb 2 = {xb ~ Xe) + (xe - Xe 1 - Xe 2 ) + {xe 1 - Xb,) + (xe 2 ~Xb 2 ), then 

\\XB - XB, ~ XbJli - \\XE ~ XE, ~ XeJli < \E\D(E) + |^i|-D(^i) + \E 2 \D(E 2 ). 

Also, note that \\xe - Xe, - Xe^Il 1 = \ E \ ~ \ E i\ ~ \ E 2\ = \B\ - \Bi\ - \B 2 \ because E 1 t~)E 2 = 
and EiL) E 2 C E. The result follows by Jensen's inequality applied to the map t \— > t q . □ 
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Lemma 3.7. Let B,Bx,B 2 be measurable subsets ofM. n . Then 

\\XB - XB 1 ~ XB 2 \\ Ll ~ {\B\ - \B X \ - \B 2 \) = 2(|5i| + \B 2 \ - \B n (B 1 U fl 2 )|) (3.15) 

= 2(15! \ B\ + \B 2 \ B\ + \B n B x n 5 2 |). (3.16) 

Proof. Consider, first, the elementary relations 

|Bi\fl| = \Bi\-\BHBi\, i = 1,2. (3.17) 
|Bn(BiUB 2 )| = |-BnBi| + \BnB 2 \ - \BDB 1 C\B 2 \. (3.18) 
\B\(BxUB 2 )\ = \B\-\Bn(BxUB 2 )\. (3-19) 

From ([3TT1) and (f37T8l) we obtain 

|5i \ J3| + \B 2 \ B\ + \B n Bi n J3 2 | = + |-B 2 | - |J3 n (Si U J3 2 )|. (3.20) 

From ([339]) and (pT20j) we obtain 

\B \ (Bx U B 2 )\ = \B\ - (\Bx\ + |J3 2 |) + (\Bx \ B\ + \B 2 \B\ + \BnB x n B 2 \). (3.21) 



Decomposing R n as LL, ai , a2 e{o,i}{ y : (xb,Xb 1 ,Xb 2 ) = (a,ax,a 2 )} we find that 

llXfl " Xfi! " XB 2 \\ Ll =\B n n ^1 + \B \ (B t U B 2 )\ 

+ 2|(Bi n B 2 ) \B\ + \(Bx\B)\ B 2 \ + \(B 2 \B)\ B x \. 

Since \(Bx n B 2 ) \ -B| can be seen either as \(B X \B)C\B 2 \ or as |(jB 2 \ B) n 

llXfl " XBi " Xftllii = \B n 5! n B 2 | + \B \ (Bx U B 2 )\ + |Bj \ S| + \B 2 \ B\. 
Using (I3T201) and (l3J9|) we obtain flUHg) ; from (pT2TT) we obtain ([335]) . □ 
From (I3.15h we see that the minimization problem in the conclusion of Lemma l3,6l is equivalent 

to 

max{|£n (Bx U B 2 )\ : B, B x , B 2 balls of radii R, Rx, R 2 }, (3.22) 
where R, Rx, R 2 are such that \E\ = uj n R n , \Ex\ = oj n R x , \E 2 \ = uj n R 2 . 

Lemma 3.8. Suppose < Rx,R 2 < R < Rx + R 2 - Then ()3.22|) admits a solution, unique up to 
isometries of the plane, characterized by the facts that: 

i) the centres of B, Bx, B 2 are aligned 

ii) ^ BxDB 2 C B, Bx(tB, and B 2 <£ B 

Hi) dB DdBx, dBxCidB 2 , and dB 2 ndB are ((n — 2) -dimensional) circles having the same radius 
(or, if n = 2, the common chords between B and Bx, Bx and B 2 , and B 2 and B all three have 
the same length, see Figure^p,). 
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Figure 7: On the left: optimal choice of B, B\ and B 2 in (I3.23p . with h = h\ = h%- On the right: 
sublevel sets \<f> < h} and {ip > h} in the proof of Lemma 13.101 (as h increases the level sets move 
along the slab S, in the direction of e). 



In addition, the solution to (|3.22j) is such that 



. 1-1 

in— 1 



on— 1 / r> r> \ -"2" 

isn^n^l^^— (.Rx + ^-i?)^ • (3.23) 

n! \R 1 + R 2 J 

The proof of Lemma 13.81 uses the auxiliary Lemmas 13.91 and 13.101 As mentioned in Section [2j 
we write a A b to denote the exterior product of a, b G M. n . In particular, we use that |a A b| = 
|b| dist(a, (b)). The purpose of Lemma 13.91 is to show that -B(p + he,R) can be written as the 
intersection of the two sets in Figure Eb), for all h G M. We then write the derivative of the area of 
the sublevel sets with respect to h as a surface integral on dB(p + he, R), using the coarea formula 
(Lemma EIUD. 

Lemma 3.9. Let R > 0, p G W l , e G S n_1 . £>e/me 



(y) := (y-p)-e-^R 2 -|(y-p) A e| 



2 



</>(y) := (y-p)-e + v / « 2 -|(y-p)Ae|2 

m £/ie infinite slab S := {y G IR n : |(y — p) A e| < R}. Then, for all h G R, 

B(p + he, R) = {y G S : 0(y) < /t} D {y G 5 : ^(y) > h}. 

Proof. By Pithagoras's theorem |y — (p + he)\ 2 = |(y — p) • e — h\ 2 + |(y — p) A e| 2 . Then 
y G i?(p + he, R) if and only if y G S and |(y — p) ■ e — h\ < \J R 2 — |(y — p) A e| 2 , that is, if and 
only if 

y G 5, (y - p) ■ e > h and 0(y) < h, 
or y G S, (y — p) ■ e < h and VKy) > h. 

This proves that B(p + /ie, R) C {<j) < h} D {ip > h}, 

{^<h}\ B(p + he, R)c{y€R n : (y - p) • e < h} 
and {ip > h} \ B{p + he, R) C {y G R n : (y - p) • e > h}. 
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From this we see that {(ft < h} n {i/j > h} C B(p + he, R), so the conclusion follows. □ 
Lemma 3.10. Let p G R n , R > 0, E C M. n measurable, and suppose that 

U n ~ x (dB{p, R) n dE) = 0. (3.24) 
Then the map y i— > |-B(y, -R) n £/| is differentiable at y = p wii/i gradient 



D y (\B(y,R)nE\) = / -yidW"" 1 ^). 

* H JdB(p,R)r\E X 

Proof. Given e G S™^ 1 arbitrary, let 0, and S be as in Lemma 13.91 By definition of <fr and ip, we 
have that 4>(y) < (y — p) • e < ?/>(y) for all y £ S, hence 

(y - p) • e < h => <f>(y) <h and (y - p) ■ e > h V(y) > & 

for all /i£K. Thus, < h} U > /i} = S and is independent of h. From the elementary relation 
\EDSinS 2 \ + \ED (Si US 2 )| = | £7 n Si | + |£nS 2 | we obtain (first for the case \EnS\ < oo, then 
for all measurable sets) 

\E n B(p + he, R)\-\En B(p, R)\ 

= (|^n < h}\ + \En{4> > h}\ - \EnS\) - (|£n{0 < o}| + \En{4> > o}| - \EnS\) 

= \En{0<<f><h}\-\En{0<i/f<h}\. 
Writing y G S as p + Ae + fie' , with |e'| = 1 and ele', a direct computation shows that 

D<p(y) = e - —=g== and £ty(y) = e + 



Hence, by the coarea formula and Pithagoras's theorem, 

\EnB {v + ke,R)\-\EnB {p ,R)\= C { f - I Z^) * 



r}nE \D4>{y)\ 7{^=r}n£; \Di/>(y)\ 
sgn(A - t) 

J dB(p+Te,R)nE 



[ [ sgn(X-r) VR H dH n ~ l (y)dr 

JO JdB(p+re,R)nE n 

= e.[ [ ^d^(y)d, 

JO J dB(p+re,R)nE U 

Since h and e are arbitrary, the above equation expresses that for all h G M n 

\E n B(p + h,R)\-\En B(p, R)\ = h • C [ ^_P XE _ rh ( z ) d^- x (z) dr. 

■/ JdB{p,R) u 

Denoting |{r G (0, 1) : z + rh G E}\ by q(z, h, E), Fubini's theorem gives 
E n B(p + h, - \E n B(p, - h • /" 5j__£ d'H n ~ 1 (z) 

JdB(p,R)r\E R 

< |h| / | Xs (z)-a(z,/ i , J E)|d^ 1 (z). 
JdB{ P ,R) 
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Due to the connexity of the line segment joining z and z + h, if dist(z, dE) > |h| then either 
z E Int E and a(z, h, E) = Xe( z ) = 1, or z E R n \ E and a(z, h, E) = Xe( z ) = 0. Therefore, 



limsup|h| 1 



EnB(p + h,R)\ - \EDB(p,R)\ -h- / ^— -dft n_1 (z) 

JdB(p,R)r\E 

< limTT-^zE dB(p,R) :dist{z,dE) < |h|}) = 'H n ~ 1 {dB{p,R) C\dE), 

h— >0 



completing the proof. □ 

Remark 1. The example p = 0, J? = 1, £ = (-1, l) n \ 5(0, 1) shows that |-B(y,#) n £| is not 
always differentiable with respect to y if (|3.24p is not satisfied. However, this condition holds in 
the situations to be considered in the sequel, namely, when E is a ball, the union of balls, or the 
intersection of balls of radii different from R. 

Proof of Lemma \3.8i The existence of solutions to (|3.22p can be easily deduced from the continuity 
of \Bn(BiL)B 2 )\ with respect to the centres of B, B\, and B 2 . Let {B , B\, B 2 ) be one such solution. 
We divide the proof of ITT) - HIT|) in the following steps. 
Step 1: one of the following possibilities occur 

dist(Bi HB 2 ,B) > 0, dist(Bi n B 2 ,R n \ B) > 0, or B 1 n B 2 = 0. (3.25) 



Suppose, looking for a contradiction, that neither B\ n B 2 f]B = nor B\C\ B 2 C -B. Then, by the 
connexity of B\ n -B2, there exists xo S Bi (1 B2 n dB. Let = i?(p, R), e := |*°Ip| , and consider 
the following parametrization of dB(p,R) using spherical coordinates 

f(0,£) : = p + (Rcos8)e+ (R sin 6)£, 0e[O,tt], £ E §™ -2 := S™" 1 n (e)^. 

Applying Lemma 13.101 to E 1 = B\ U -B2 (see Remark [I]) 



-£-(|B( P + fce,i2)n(BiUB 2 )|) 

dn, 



z ~ P 

h=0 JdBn(BiUB 2 ) R 



e 



= R n ~ 1 / cos0(sin0) n - 2 d^dH n - 2 (O 

J§L l ~ 2 J0£(0,n):f(0,£)£E 

We can write the integral with respect to 8 as 

cos #(sin0r- 2 ( X £(f (#,£)) " XiKffr " 

JO 

If we prove that 

f(vr - 0,£) E T^uTh f(0,£) E STuBa" for every 9 E [0, tt/2] (3.26) 

and that 

X.e(f (#,£)) — xe(f( 7r ~~ ^)^)) = 1 f° r an (0>£) m a se t OI " positive measure, (3.27) 
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we will obtain that (|£(p + he, R) n (£1 U B 2 )\) > at h = 0. The contradiction will follow by 
noting that if (B, £1, £ 2 ) solves ([332]) . then £ x |£(x, £) n (£i U £ 2 )| must be zero at x = p. 

Suppose that {(n — 9o,£) G Bi for some i = 1,2 and some #o G [0, Since £n<9£ is connected 
and contains /(0,£) = xo, its projection to the plane p + (e, £) must contain the whole of the arc 
f (#,£)> ^ G [0,vr - 6> ). This proves (I33HD . In order to prove flQZP , define 6>i(£) := sup{6> G [0,vr] : 
f(0,£) G £i U £ 2 }. Arguing as before, we see that 

\{e€[0,*]: X E(f(e,Z))-XE(f(n-e,Z)) = l}\>0 (3.28) 

unless 0i(£) = or 6>i(£) = vr (by continuity, if (EQgJ) holds for at least one £ G §™~ 2 , then (13371) 
follows). Since R\,R 2 < R, in fact 6\ = ir is not possible (in that case xo and xo — 2Re would 
belong to some Bi, but diam£ = 2Ri < 2R). It remains to rule out the possibility that = 

for all £, that is, that B n (£i UB2) = {xo}. If that were the case then B and £1 would be tangent, 
so for all h < R\ we would have that 

|£(p + he, R) n (£1 U £ 2 )| > \B{p + he, R) n Si| > = \B n (£1 U B 3 )| 

and {B,B\,B 2 ) would not be a solution to (|3.22p . This completes the proof. 

Step 2: the centres of B, B\, B 2 lie on a same line. In all the three cases considered in (|3.25p . 
\B n Bi n B 2 \ = \ (B + h) n B\ n B 2 \ for every h sufficiently small. Also, for given R, R 1 , R 2 , the 
expression \B(yi, R4) D B(y, R)\ is a decreasing function of |y — yj|, i = 1, 2. If y were not in the 
line containing yi and y 2 , both |y — yi| and |y — y 2 \ could be reduced by displacing y towards 
that line. By (|3.18|) . this would increase \B Ci (Bi U B 2 )\, contradicting the choice of {B, B\, B 2 ) as 
a solution to (13.221) . 

Step 3: {B,B\,B 2 ) satisfies ^ '-[m p . Moreover, these conditions uniquely determine the distances 
and relative positions between the centres (that is, the solution to (|3.22p is unique up to isometrics). 

Let h, hi, and h 2 denote, respectively, the radii of dB\ n dB 2 , dB n dB\, and n 9i?2 (or the 
semi- lengths of the common chords between B\ and B 2 , B and B\, and B and £2 if H = 2) defining 
these radii (or lengths) as zero in case of empty intersection. By virtue of [1]) , both pi — p and P2 — p 
are parallel to e := ipgZpj ; where p, pi, P2 are the centres of B, B\, B 2 , respectively. Setting 
qi := (pi — p) • e, i = 1,2, and using Cartesian coordinates (yi, . . . ,y n ) with p as the origin and 
e in the direction of the yi-axis, we have that B = -B ((0, 0, . . . , 0), R) , B\ = B((qi,0, . . . , 0), i?i) , 
B 2 = B((q 2 ,0, ... ,0),R 2 ). By (I3T8D ancfl Lemma EM 

is n (Si u s 2 )| = A|B n B t \ - ^-\(b n s 2 ) n B x \ 

oqi oqi dqi 

Zx ~ qx me-\ zi ,...,z n )- [ me-\z 1 ,..., Zn ). 



— v . i5 , „„, . 
dBillB u i J9_Bin(Bns 2 ) n i 

In the first of the possibilities considered in (|3.25p . B cannot intersect both B\ and £2, hence 
(B, Si, £2) is not optimal (for example, it would be better if B contained completely either B± or 



8 There is exactly one situation not covered by Lemma 13.101 namely when Ri — R2 and B\ = B2 CC B, but it is 
easy to see that this does not give a maximum of \B PI (Bi U i?2)|- 
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B2). In the other two cases we have dB\ n (B n B2) = dB\ n -E?2- Parametrize <9I?i by 

ze^ O z- Pl = cos 0)e + ( J R 1 sin#)£, 0€[O,tt], £ G §™~ 2 := S"" 1 n (e) x . 

By definition of e, q\ < q2- Therefore, z G (9i?i n 2?2 if and only if 8 G [0, 6*2), where #2 is one of the 
two angles in [0, it] such that by h = R\ sin 62 (when h = 0, we choose 62 = or 02 = vr according 
to whether 5 2 (1 Bi = or B 2 C Bi). Thus, 

^-((Bn^nfiil =^ n ~ 2 (§r 2 ) r 2 J R n ^ 1 cos0(sinfl) n - 2 d# = w n _i/ l n - 1 . 

As for the integral on dB\ n -B, the same argument shows that it equals — (sgngi)a; n _i/j" _1 . Af- 
ter obtaining the corresponding expression for -£^\B D -B2I, and by virtue of the optimality of 
(B, Bi, B2), we obtain 

sgn^)/^- 1 + h n ~ x = h n ~ l - sgn(g 2 )/i2~ 1 = °- 

The case h = h\ = Ji2 = is not optimal (due to the assumption R < R\ + R2), hence qi < 
< q2 and h = h\ = /12 > 0. This proves ln|)-liii|). It remains to show that qi, ^ and h are 
uniquely determined by these conditions. Denoting the hyperplane containing the intersection 
of the boundaries of two (intersecting) balls B',B" by H(B' , B"), we have that the hyperplanes 
H(Bi,B), Tl(Bi, B2), and H(B2,B) are given by {yi = a±}, {yi = a}, and {y\ = 02}, for some a±, 
a, ci2 G R. Clearly, the following must be satisfied 

(01 -qi) 2 + h 2 = R\ (a - qif + h 2 = R\ a\ + h 2 = R 2 

a 2 + h 2 = R 2 (a - q 2 ) 2 + h 2 = R\ (a 2 - q 2 ) 2 + h 2 = R 2 2 . 

In particular, |ai| = \a,2\ = V R 2 — h 2 , \a\ — q\\ = \a — qi\ = \J R\ — h 2 , and \a — #2! = \ a 2 — Q2\ = 
\JR\ — h 2 . Conditions lu])-lm]) imply that a\ < q± < a < q2 < (12 and a\ < < 02- Therefore 

qi = ^R 2 -h 2 - VR 2 ~ h 2 , q 2 = VR 2 ~h 2 - yjRl - h 2 , (3.29) 
which shows that q\ and (72 are determined by h. We also find that 



a-q 1 = ^R 2 -h 2 , q2-a=^R\-h 2 . (3.30) 
Adding the equations in (I3.30P and subtracting the equations in (|3.29p yields (see Figure [8]) 



q 2 -qi = VR 2 ~ h 2 = ^R 2 -h 2 + y/R 2 - h 2 . (3.31) 
We may assume, without loss of generality, that R2 < R\. Rewrite (|3.3ip as 

R 2 -R 2 



R 2 -h 2 = Q. 

VR 2 -h 2 + s/Rj - h 2 

The expression at the left-hand side is increasing in h, and equals R — (Ri + R2) < at h = 0, 

R 2 — R 2 

and , — j > at /i = i?2- This shows that h is uniquely determined by i?, i?i,i?2, and 

hence the balls B\ , B2 too. 
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Q2 - gi I ^ 1 

1 — — 1 2 ^ 1? 

Figure 8: Relationship between h and the distance between the centres. 

Step 4 : proof of (|3.23j) . For each k G {2, . . . ,n} denote by the /c-dimensional polyhedron 
with vertices (the convex hull of) 

{(q 2 - R 2 )e, (qi + Ri)e} U {ae ± he, : i = 2, . . . , k}, e; := (0, . . . , 1, . . . , 0) . 

V * ' 

j-th position 

It is easy to see that T-L 2 {P 2 ) = hj, where 7 := \(q\ + R\) — (q 2 — ^2) I , and that 1~L k {Pk) = 
2KH k - 1 {P k ^i)/k, for k € {3, . . . , n}. Thus, \P n \ = (2 ri - 1 /i"- 1 7 )/n!. 

From the previous analysis, we have that B\ n B 2 contains P n . From this we obtain (|3.23p . 
since, by virtue of (|3.3ip . 



7 = i?i + R 2 - VR 2 ~ h 2 > Ri + R 2 - R, (3.32) 

and 7 = £_ + ^_< ( ^ + ^ 2 . (3.33) 

i?i + V^i - ^ 2 R 2 + y/Rj^f? R1R2 y ' 

□ 

We finally prove the main result. 

Proof of Proposition li.ffl We can assume that + I-E2I™ > (otherwise the estimate is 

trivially true). By (|3.16p and (|3.23p we have that 

n-l 

iyn / d td \ 2 

min(||xB - Xfli " XB 9 |Ui " (|5| - |-Bi| - |£ 2 |)) > -r(i?i + ife - fl) V 

n! \Ri + R 2 J 

where the minimum is taken over all balls B, B\, B 2 with \B\ = \E\, \B\\ = \Ei\, \B 2 \ = \E 2 \, and 
R,Ri,R 2 are such that \E\ = uj n R n , \Ei\ = 0J n R±, \E 2 \ = uj n RV,. Thus, by Lemma I3T6I 

n n n n+1 n " 

\E\D(E)~ + \E 1 \D(E 1 )~ + \E 2 \D(E 2 )~ (R x + R 2 - R)~~ ( R X R 2 



\E\ + \E 1 UE 2 \ (R n + ff{ + R%)^ \Ri + R 2 

The quantities R n + i?" + RV^ , R\ + R 2 , and + /?2) n are comparable, since we are assuming 
that R < Ri + R 2 and by virtue of the identity a n + b n < (a + b) n < 2 n ~ 1 (a n + 6 n ). Hence 

n rP 1 n(n-j-l) n 

(iT + .R? + i?£)— < (7(^1 + ^2)^ =C(i2i + i? 2 ) 2( - 1) (^1 + ^2)2, 
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which implies that 

n n n n(n+l) n n 

\E\D{E)— + \E 1 \D{E 1 )— + \E 2 \D(E 2 )— > / Ry + R 2 - R \ R'{ R% 

\E\ + \E 1 l}E 2 \ ~ V R1+R2 J (Ri + R 2 ) n ' { ' ' 

By the mean value theorem, there exists £ between R and R\ + R 2 such that 

f? _l r _ d _ ( Rl + R ^ n - Ri - R 2 f (Ri + R2T - -R" 

Since we are assuming that R < R\ + R 2 , then £ < i?i + R 2 and 

Ri + R^-R > 1 {Ri + R^j-^-^ f (lE^ + \E 2 \^) n -\E\ \ 
R1+R2 ~ n {Ri + R 2 ) n \{\E l \^ + \E 2 \^) n - \E X U E 2 \) ' 

Suppose now that |£i| > 1^2 1, so that R ^ R2 > ^- By the binomial theorem, 

{Rl+ R 2) n_ R n_ R n = ^ / X / \ re-fc / \* _R^_ 

(#1 + ^2)™ ^WU + W \Ri + R2j ~2 n ~ 1 R 1 +R 2 l " j 

(we have considered only the term corresponding to k = 1). Combining (|3.35|) with (|3.36|) we obtain 

n(n+l) n(n + l) / 1 1 \ + 

R 1 + R 2 -R\*^ ( R 2 [ i\E x \- + \E 2 \~) n -\E\ \ 2(n_1) 



-Rl+#2 / \R 1 +R 2 J \{\E X \n + \E 2 \n) n - \Eil] E 2 

Ri 

R1+R2 — 2' 



The conclusion follows from ()3.34p and the above equation, considering that B > I - D 



4 Upper bounds 

As explained in the Introduction, we obtain the upper bounds of Theorem [2] and Corollary Q] by 
finding suitable test functions opening cavities of different shapes and sizes, the main difficulties 
being to satisfy the incompressibility constraint and the Dirichlet condition at the boundary. We 
split the problem into two: in Section 14.11 we define a family of incompressible, angle-preserving 
maps whose energy has the right singular behaviour as e — > 0, with leading order {y\ + v 2 ) \ log e| , 
and serves to define the test maps close to the singularities. In Section [4.21 we extend those maps, 
using the existence results of Riviere & Ye [63]) in order to match the boundary conditions. 

4.1 Proof of Theorem [2] 

In order to compute the energy of the test functions, we will need the following auxiliary lemmas, 
whose proof is postponed to Section 14.31 

Lemma 4.1. Let VL be a domain in W 1 , star-shaped with respect to a point a € W 1 , with Lipschitz 
boundary parametrized by £ 1— > a + q(C)C> C £ § n_1 - Let v > and define u : R n \ {a} — > W 1 by 

u(a + rC):=Aa + /(r,C)C, f(r, CT := r n + (A™ - l)q(C) n , r€(0,oo), C £ § n ~\ (4.1) 
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with A n := 1 + tjjt. TTien u is a Lipschitz homeomorphism, detDu = 1, u(x) = Ax for all x G dfl, 
u(H\ {a}) = AO\im T (u,a), u(M n \ 0) = K n \ AO, |im T (u,a)| = u, and /or aZZ r, £ 

ra-1 



„n-l 



Du(a + rC) 



n 



1 



<C r + H» 



max{q, |£>q|} \ 



+ 



q(C) n maxjg, iDgl}"- 1 !^^ v 

|n| |n| y r : 



C freing a constant depending only on n. 

Lemma 4.2. Suppose that a G R n , < d < p, and a = a+de /or some e G § n_1 . Zei £ i-> a+q(£)£, 
£ G S n_1 be the polar parametrization of 8B(sl, p) taking a as the origin. Then 

2 



i) for all C G § , |g(C)| < 2p, |D?(C)| < 2d|C A e|, and |£>g(C)l < 2d 



g(C) 



m) i/ C • (a — a) < i/ien g(C) > p |C • e| and 1 < 



?(C) 



d|C-e| + ^(/o-d) 



Vp(p~ d ) 

< 2 



|C A e| 



z/ C • (a — a) > then 



< 



2^2 



VpJp - d) 1 + 



Lemma 4.3. Let < d < p, a G R n , e G § n ~\ and := {x G £(a,p) : (x-a)-e > p-2d}. T/ien 

n|Sl| > uJn-id 2 (2p — d) 2 . 
Proof of Theorem^ - Step 1: Construction of the domain. 

Let ai, a2 G R n and d := |a2 — ai| > 0, as in the statement of the theorem. Call e := r~Efj • 
Given di, d2, pi, and p2 such that < d\ < p\ < d, < d2 < P2 < d, and d\ + d2 = d, define 

ai := ai + (pi - di)e, a 2 := a 2 - (p 2 - d 2 )e, Bi := J3(ai,pi), B 2 := B(& 2 , p 2 ) 

(ai, a2 are chosen such that BiL)B 2 fits in an infinite slab of width 2d, as in Figure [S]). As stated in 
the Introduction, our aim is to show that for every 5 G [0, 1] there are unique di, d 2 , pi, and P2 such 
that the ratio between the width of B\ n B 2 and that of B\ U B 2 is exactly 5 (i.e., <5 := 



and such that JS4 = — , with 



d 



fil := (x G .Bi : (x — ai) ■ e < di}, 2 := {x G -B 2 : (x - a 2 ) • e > — d2}. 



(4.2) 



To this end, we will first consider a simplified but equivalent problem. Fix d > and e G S n , 
and let 5 := {x G W 1 : |x • e| < d}. Given pi and P2 in (0, d) define 



B i = B({-d + pi)e,pi) and B 2 = B((d - p 2 )e, p 2 ) 



(4.3) 



(the balls of radii pi, P2 contained in S and tangent to dS from the right and from the left). If the 
balls intersect, let a G (— d, d) be such that x • e = d for x G B\ n B 2 and define 



rii := {x G B\ : x • e < a}, Q 2 := {x G -B2 : x • e > a}, 
We want to show that 



Pmir 



-ufd 



(4.4) 
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|Oa| 



ii) for every p\ G [p m i n , d) there exists a unique p2 G [0, d] such that i?i n B2 / and -U^ 



C-2 
Vl 



iii) p2 = p2(pi) is such that P2 < pi and such that the ratio Pl+ ^ 2 rf increases from to 1 as p\ 



d 

increases from p m i n to d. 



This will imply that for every 5 € [0, 1] there are unique p\ and P2 such that Pl = 5 and 

= ^. Let d\ := and d 2 := with d as in (|4.4p (they are the semi-distances from the 
plane containing B\ n i?2 to the walls of the slab S containing B\ U -E>2)- Based on the previous 
reasoning, it can be seen that these values of d\, d 2 , pi, and p2 constitute a solution to the original 
problem, and that no other choice is possible. 

In order to prove 0), define B[ := ((-d + p min )e, p min ) , B' 2 := (p min e,d- p min ) {p min is such 



that B[ and B' 2 are tangent and = ^). If < p x < p min and B 1 n B 2 / 0, then pfj > ^, 

since fii C -B^ and VL2 D -B 2 - Hence -^j = ^ => pi > p m i n) as claimed. 

Fix pi G [pram,d). In order for i?2 to intersect B\ we must have that P2 > d — p\. When 



P2 = d — pi, fii and are tangent balls with IS 2 - = l " < l " ^ mm,) = ? < 1- It is clear that 
| fill decreases and 1 increases as p2 increases (the intersection plane moves to the left), therefore 
j^j is increasing in p2. When p2 = pi, the ratio is 1. This proves [n]) and the first part of|m|). A 
similar argument shows that pl+ P^~ d [ s increasing in pi (it follows from the fact that if we fix p2 
and increase pi then the intersection plane moves to the right and decreases). 

It is clear that if pi = p m ; n then P2 = d — p m ; n and P 1+ P 2 ~ d = 0. It only remains to prove that 
as pi — > d also p 2 — > d. By (|4.3|) . \B 2 \ B x \ < \B(0, d)AB x \ ->■ as pi d, hence 

lim W = lim I" 1 ' (l + mi^MlM ] = v .± (i i ^Pi^d\B 2 \B 1 \ \ _ Vl 



For pi < d, dBi n dB 2 is of the form A(pi) := {a(p 1 )e + yj p{ - a(pi)V : e' G S n_1 ,e' _L e}. 
Since a(pi) is determined by j^jj, it has a well-defined limit as pi — > d. The sphere <9i?2 can be 
characterized as the one containing A{p\) and the point de. Thus, in the limit, dB 2 will be the 
sphere containing de and A(d), which is none other than dB(0, d). In particular, P2 — > d, as desired. 

- Step 2: Definition of the map. 
We define u : W 1 \ {ai,a2} piecewise, based on Lemma 14.1} in the following manner. Inside Vt\ we 
apply Lemma 14.11 to fi = Vt\ and a = ai; inside £l 2 we apply Lemma 14.11 to Q = Q 2 and a = a2- 
Finally, in order to define u in M n \ Cli U £l 2 we define 

(ai + pie) + (a 2 - p 2 e) _ 
a = = ai + (d - p 2 )e = a 2 - {d - pi)e 



(when 5 = 0, a* is the intersection point; when <5 = 1, a* is the center of the ball) and use Lemma 
Owith ft = Sli U n 2 , a = a*. Let C ^ »i + <?i(C)C, C ^ a 2 + <? 2 (CK, and C ^ a* + q(C)C be, 
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respectively, the polar parametrizations of dfli, 80,2, an d d(Q\ U f2 2 ) (with ( G §" 1 in all cases). 
To be precise, 



u x 



Aai + |x-ai| n + 



Aa 2 + (|x-a 2 |" + ^g 2 



x— ai \ \ " x— ai 

l x ~ a i| / J i l x - a ll 

X— SL2 \ n \ n X— SL2 

|x-a2| J J |x-a2| 

Aa* + ( |x-a*r + - 5 



x G J2i \{ai} 
x G H 2 \ {a 2 } 



with 



A n -1 



V 1 + u 2 



'■'2 



|f2i| \Q 2 \ |f2iUfi 2 | 

Since 1^4 = the construction is well defined and u(x) = Ax for all x G 8fl\ U9f2 2 . The resulting 
map is an incompressible homeomorphism, creates cavities at the desired locations with the desired 
volumes and is smooth except across <9f2i U <9f2 2 (where it is still continuous). It only remains to 
estimate its elastic energy. 



- Step 3 : Evaluation of the energy in W 1 \ (fix U ^2)- 
By Lemma l4.2ip . max{q, \ Dq\} < 2d, then, by Lemma |4. II 



-n-l 



Du(rQ 



n <C r+ d{Vl+V2)t 



1 \ n—l 



If2i U f2 2 



+ 



+ 



Cd 



n— 1 



fil U f2 2 | 1^1 U f2 2 | / r 



5) 



1 



Since pi, i = 1,2 increases with 5 and assumes the value — r i — j- when 5 = 0, it follows that 



2oj n d n > u n (p? + fi) > U Q 2 \ > ^^(Pl + P2) > 2~ n u; n( r 



(4.6) 



(since f2iUf2 2 D 5j for each i = 1, 2). Consequently, for any i? > (using that / H n (§> n ) = nu n ) 



1 

n -/B(a*,i?)\QiUf2 2 

~T~ X C 

< 



Du(rC) 



1 



dx= - 



max{g,i?} 



„n-l 



(w nJ R n + 2>x + W2 )) + ( Wl + t; 2 ) 



L>u(rC) 



drd^-^C) 



R 



dU 



n-l 



where (logx)+ := max{0, log x}. Note that ^log^J < (log^) + + [log -J . Also, 

|Oi U Q 2 | = /* f 9C r n - 1 drdn n - 1 (C) = -f oo n q{Q n d% n ~ l {C). (4.7) 
Finally, since |a* - fii| + |a* — a 2 | = d(l — 5), Lemma [^73T|) implies that \Dq\ < 2d(l — 5). Hence, 



n 



B(a*,R)\QiUn 2 



Du 



dx <C(vi +v 2 + oo n R n ) + ( Vl + v 2 ) (l + C(l - 5)) ^log ^ 
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Figure 9: Angle 8q > ^ and choice of spherical coordinates for 5 = 0. 



The main problems at this point are that if 5 — > then p2 is of the order of 



x 



X ISO - 



on dB2 n 80,2 if ^ — > 0) and q(£) tends to vanish on dB\ n <9i?2 (see Figured]). Parametrize S 
by C = -cos(9e + sin6>£ with 8 G (0,vr) and £ G S* := S"^ 1 n (e) 
only study the term with \Dq\, that is, we are to prove that 



oo 

n-l 



Since log -| is bounded we 



H 



n-2 




sin#) r 



. 2 IAz(C(fU)l 



d 



log 



d 



dO 



is bounded independently of d, 5, v±, and V2- It can be shown that a* + q(9, £) £ <9I?i for all 

9 G (0, 5) (due to the fact that pi > p 2 , see Figure [9]) , and clearly £ ■ (a* — Si) = — cos 9 (d — p2 ) < . 
Lemma l4.2iip can thus be used to estimate the first integral by 



2 ft, d 

— log ■ 

a pi cos I 



d9 <2[ max \t log £ 

\te[o,i] 



log 



1 /7T 



2 V2 



d# 



1 + log 



7T 



As for the second integral we divide (f , 7r) into (^, #o]U[#o> t 1 ")) according to whether a*+q(9, £)C(# ; £) 
belongs to <9i?i or to dE>2- For 9 > 9q we can still use Lemma l4.2tip (this time with a = a 2 and 
P = P2) to obtain exactly the same upper bound as before. For 9 G (|, #0), use parts []} and[m|) of 
Lemma 14.21 together with p\ — |a* — 3l\\ = d5 to obtain 

-2 

sin#. 



2{d- P 2)q 2 . n , , .A (d - p 2 )| cos 01 

1 < v _ r / -^sine and |J3g < 16(d - p 2 ) 1 + 



d 



<5pi d 2 



\JJpx~d 



Then, for any a G (0, |), using that £ 2q | log i | < (2ae) 1 for every £ G (0, 1), 



0o 



IAzI 



log- 



d0 < 



r0 


Dq 


1-0 




K) 


k 


d 




d 



< 



< 



2 a 16 



l-a 



2ae 

2 l-a /fy 



d - P2 



l-a 



d0 

d - P2 
Spi 



1 + 



(d - p 2 )| cos I 
y/Spid 



2(a-l) 



sin#d0 



ftp 



d 



d~P2 



/Sp-^d 



1 (l + t) 



2a-2 



d£ 
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The last integral can be bounded by means of the relation 

1 



(1 -2a) [ X {1 + t) 2a ~ 2 dt = 1 
J o 



< 1 



1 



(l + x) 1 - 20 ~ l + x 1 + x' 
Using that 7 + ^fT^j > 1 for all 7 G (0, 1) (applied to 7 - ^ - |a *~~ ai1 ' 



pi 



pi 



ll—a 



log - d9< 



a(l — 2a)e \ d 



Spi\ 2 a d - p2 



1 



Pi 7 + V 1 -7 



)1 — Q 



< 



< 



a(l — 2a)e 

gl-a 

a(l — 2a)e 



52" 



d- P2\ 2 a ( d- p 2 



d 



in 



5i- a (l-8)v 



We conclude that for all R > 



1 

n 



B(a*,R)\n 1 un 2 



dx < C{v 1 +V 2 +OJnR n ) + {vi+V 2 )[l + C{l- 5)) ( log — 1 . 



- S'iep ^: Estimating the energy in 
Near the cavitation points we still have that f oj n qf d7i n = |f2j|, i = 1,2, so by Lemma 14.11 



1 

" Jcii\B e .(in) 



Du 



dx < Civi + w n/ o™) + «j log 



2 Pi 



Si 



+ C 



Vi + V 2 
U Q 2 



f maxfe l^l}™" 1 !^! d^ 1 ) log 
Js™- 1 J 



2d 

Si 



< C( Vi + LOnpf) + Uj log — + C( V1 + V 2 )^ 

Si. d 



riff \m\ log ™ 



For Cli set C = — cos#e + sin 6>£. If 6 E (0, f ) then, by Lemma [472], using that |ai — fii| = pi — d\, 



\D qi \sin n - 2 6de< 16(pi - di) / " (1 + 



Pi - di 



y/dipi 



cos 6 sin#d# 



pi- d i 
V c1 ipi 



(l + t)- 2 dt 



di pi-di 
pi7 + V / l- r 7 : 



with 7 = 1 - & Since 7 + v / T Zr 7 < 1 for 7 G [0, 1], 



^r 1 



\D qi \ S m n - 2 ede < P i- 2 y/d7ft{pi -d 1 ). 



Define 9\ as in Figure [5j By Lemma 14.21 |-^<7i| ^ 2(pi — d\)sm.9 and (/1 > \/dip\, hence 
Pl [ 1 |Dg 1 |sin n - 2 0d0< pi(/9i-di)| cos 6> x | < { Pl -dx)^- < y/d^(pi - di). 
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For 9 € it), qi(C) is given by qiC • e = di hence 



Qi(0) 



di 



COs(7T — 6) 



and \Dqi(C(9^))\ 



gi(i-C®C)e 



-C-e 



dl sin 
COS 2 (7T — 9) 



Using that 1 — | cos 9i \ 



sin 2 gi 
1+| cos 01 1 



< sin 2 9i and that > (pi — di) cos + \fd~ip~i > y/dip~i, 



Pi / |£>gi|d0<di/9i 



di di sin" 5 0i 

cos^l * 2 ~ Pl cos(7r-6»i) 



pi(gi(fli)sinfli) 2 /- 

^ <Wdipi(pi-di), 



the last equality being due to the fact that q(9i) cos 9i = di and &i+q(9i)£(9i, £) 6 dB(&i, pi). Now 
we show that max{gi, |-Dgi|} < 8pi. The fact that q(9i) > y/dipi implies that pi\ cos#i| < V^iPi- 
Clearly is decreasing, therefore 



q(0) < q(9i) < 2((pi - di)\ cos0i| + y/d^) < A^/d^i < A Pl . 



As for we have that qi(9)sm.9 is decreasing and g(#i)sin#i = 2-sJ di(pi — di), then 

\Dqi 



qi(qism9) < 2 gl (6> 1 ) > /d 1 (p 1 - d Q ^ _ ^ ^ ^ 



gi cos(7r — C/J di 
The study of u in $7 2 being completely analogous, the conclusion is that for all R > 



B(a«,R)\(B £1 (ai)UB £2 (a 2 )) 



V?Td 



R R 

dx < C(«i + v 2 + oo n R n ) + vi log h«2 log — 

£1 £2 



+ + „ 2) ( (1 - «, (,„ g f) + + ^^io g | + io, I) 

In the case of ai it is pi — di that has an interesting behaviour, whereas for a2 it is d 2 - This follows 
from our final ingredient: the 'height' of B(&i, pi)C\ B(a 2 , p 2 ), whether we measure it from the first 
ball or from the second, is the same. The corresponding expression is di(pi — di) = d 2 (p 2 — d 2 ). 
As a consequence, 



Pi ~ di 



5(pi - di) 



5do 



s d 4. 



d (pi - di) + (p 2 - d 2 ) di + d 2 d 

The theorem is thus proved since, by Lemma 14.31 



n+l 



\n 2 



n + l 

P 2 d 2 



< C- 



t' 2 |^lU^2| 



n-1 
2 



< c 



n + l 



v 2 



Vl + v 2 



-T^d d- 



□ 
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4.2 Transition to radial symmetry 

Our proof of Theorem [3] is based on the following result (see [531 EOJ E21 EU HU [8] f° r related work): 

Proposition 4.4 (Riviere- Ye, [M], Thm. 8). Let D be a smooth domain, k = 0, 1, 2, . . . and suppose 
that g E C k,1 (D) = W k+1,oc (D) with info g > and f D g = 1. T/ien, i/iere exists a diffeomorphism 
(f) : D —> D, satisfying detDd) = g in D and d> = id on dD, such that, for any a < 1, 4> is in 
C k+1 ' a (D) and 

||0 - id|| cfe+1 , Q(S) < C||g - l|| cfc ,i(D) 
/or any < 5 < 1, where C depends only on a, k, D, ini^g, 5, and \\g\\os. 



Lemma 4.5. Let £ € S n i— >• a* + q(C)C be the polar parametrization of d(fl\ U Q 2 ) define 

p(C) n ■■= R t+^ +v ^ J^l 2 y Ce^ 1 , (4.8) 

i?i being fixed and such that £li U £l 2 C B(a*,Ri). Suppose that u is a one-to-one incompressible 
map from {R\ < |x — a*| < -R2} 0^0 {?*C '■ p(C) < r < R3}, for some R2, R3 > 0. Then 

UJn (R^-R^>l^i-(v 1 +v 2 )(l-5). 
Proof. Denote maxgn-i q = 2p\ — 5d by g max - By incompressibility (using that Ti n ~ 1 (S n ~ 1 ) = nw„), 
u n iq - uj n Rl = |{x : R 1 < |x - a* I < R 2 }\ = IK : p(C) < r < R 3 }\ (4.9) 

= [ r r n ~' dr dte-HC) = ^3 - ~ fri + v 2 ) fs "-;" ^oT" 1 ■ 

Hence, the requirement that -R3 > p(£) for all £ £ S rt_1 is equivalent to 

, w « i§n-i(<7max ~ 9™) d"H n_1 



w n (^' - 1%) > ( Vl + v 2 ) 
Write C := - cos 9e + sin#£ with € [0,tt], £ € £ := S"" 1 n (e)^. For all 6 G (0, f ] 



|^i U Q 2 \ 

Pnr oil # C CO _ 



<?max — = 2pi — <5g£ — (pi — Sd) cos — \J Sd(2p\ — 5d) + (pi — Sd) 2 cos 2 9 
(pi + (Pi - Sd)(l - cos9)) 2 - (5d{2 Pl - Sd) + (pi - Sd) 2 cos 2 9) 

pi + {pi - Sd)(l - cos 9) + y/5d(2 Pl - Sd) + (pi - Sd) 2 cos 2 
(pi - o"d) 2 (sin 2 fl + (1 - cosfl) 2 ) + 2p x (pi - Sd)(l - cos 9) 
> (2pi - 5d) + (2pi - Sd) + p!-5d 

2( Pl -Sd)(2 Pl -Sd)(l-cos9) ^2^ flK 2d 
= — _ > -(d- p 2 )(l - cos<9) > y(l - 5)(1 - cos<9), 

where we have used that p\ — db = d — p 2 and P2 < d. Therefore, 

Q0nf Sn ->(q^-q n )dH^ n n-2 {s) fl( qm&x - q ) q ^( S m9r- 2 d9 

|OiUn 2 | «w n 2d" 2^ 2 3vr l h 1 j 

□ 
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Proof of Theorem [J We prove the theorem in the following stronger version (see the remark after 
the proof of Corollary [I]) : " Let R\, i? 2 be such that 

R x >2d and io n (R% - R 7 ?) > 4 n n( Vl + v 2 ){l - 5) (4.11) 

(5, v\, V2, sli, a.2, d, e±, and £2 being as in the original statement). Then there exists a* ; C\, C2, 
and u : W 1 \ {ai, — > W 1 such that u|]Rn\.B(a* t R 2 ) ^ s radially symmetric and for all R> R\ 



1 

n ys(a*,B)\(B ei (a 1 )US e2 (a 2 )) 



Du 



R R 

dx < C\(vi + v 2 + oj n R n ) + v\ log h V2 log — 

£1 £2 



+ C 2 ( Vl + v 2 ) ((I - 6) log ?± + 5 ( p/^log - + 2 f/^log - 

+ (, 1+ , 2 + ^ 2 )£^ — j ij, 

i/ie function S 6eina sac/i i/iai £(t) < 00 for t £ [0, ^) and E(i) = 0(t n(n_1) ) as t -)■ 0". The 
Theorem follows by choosing R\ and i? 2 as in (jl . 15 j) . 

Since the constant in Proposition 14.41 depends on the reference domain, we work on the annulus 
L>:={zeIR n :l<|z|< \/2} (we choose \/2 so that \D\ = u„). Our strategy is to define 
u : £?(a*,i?i) \ {ai,a 2 } — >• W 1 as in Theorem [2] and to look for a map 

u : {x £ R™ : i?i < |x — a*| < i? 2 } — > {y = Aa* + r£ : p(C) < r < R 3 , C € § n_1 } 

(where p is defined in (|4.8j) ) of the form u = v o </> _1 o w _1 , with : D — >■ D a diffeomorphism and 

w(rC) := a* + ((2 - r n )R? + (r» - 1)1$)* C, 

v(rC) := Aa* + ((2 - r>(C) n + (r n - 1)1$) » C- 

The maps w and v are parametrizations of the reference and target domains, and are defined so 
that detDw is constant and v o w _1 sends dB(a* , R), R\ < R < R2 onto a curve enclosing a 
volume of exactly v\ + w 2 + oj n R n (as can be seen by writing 

— '<*' + -0 - - + U + (l + f^if^P) ) * C (-3) 



and by considering that |{Aa* + r£ : C G S n_1 , < r < p{Q}\ = f sn _! w n p n d^™" 1 ). The problem 

det Dv 
dot flw 



for d> is d> = id on dD, det D(f> = g := ^ ^ in D. To use Proposition 14.41 we need to bound 



for all C G § n_1 , r E [1, \/2] (the constant in Proposition 14.41 depends on ||g||o,<5, so it is not sufficient 
to control only \\g — l||£°o). Using f)4.6|) and the fact that p\{S) < d and q(C) > 5d for all 5, C, 

■*,&■-,(&. - - o d«-' ge - jg - « _ 

OiUfi 2 | ~ 2-"d™ - v ; v ; 
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By Lemma SHU), 

This and Proposition 14.41 imply the existence of a (piecewise smooth) solution such that 

' ( Vl + v 2 )(l - 5)' 



\<t>- id llcM^) - S 



(4.16) 



for some function S satisfying < oo for t £ [0, ^^) and £(i) = O(i) as t — > 0. 

Define u = v o o w. Writing x = w(0(z)) and using (|4.12p and det D<p = g we obtain 



Dv(z) adj L>0(z)Dw- 1 (x) 



det£)0(z) 



i2? ll-||(?-l|| L o 



(4.17) 



Combining (|4.9p and (|4.7|) we obtain that for all R < R 2 

R 



B(a*,fl)\B(a*,Hi) ^1 



3 dx < (Ul + U 2 + Wnfij) 
(m+V2)(l-^) 



io n R n - u n R^ 



uj n Ri 



(4.18) 



By (|4TT4l) and (ITO]) . < 4"nt, with i := ZT(rI-r{) • Hence > b y ®3B, (Effl) > and ([TO]) . 

i D u ( x )r dx < ota + „ + (§ - 1 

S(a*, J R)\S(a*, J Ri) V it x J / V^l 

for R!<R< R 2 , where E(t) := gSSff » * e t ' i^) and £(*) = O^"" 1 )) as i -> 0. 

The map u can be extended to M n \ B(a* , R 2 ) by u(r£) := a* + (r n + vi + v 2 )"C- It satisfies 



1 



n JB{a*,R)\B(a*,R 2 ) 



Du 



dx < C(yi + v 2 + uj n R n ) + (vi + v 2 ) log 



R_ 

To 



The energy inside B(a*,Ri) has been estimated in Theorem El This completes the proof. 



□ 



Remark 2. For Dirichlet boundary conditions that are not necessarily radially symmetric, the above 
method can still be used provided there is an initial diffeomorphism v, from the reference domain 
D = {z : 1 < |z| < a/2} onto the desired target domain, for which g :— i)w * s bounded away 
from zero. The energy estimate will depend on info g, HDvUooHDw -1 !!^, and ||g||oo + H-^fflloo- 



4.3 Proof of the preliminary lemmas 

In this Section, we give the proofs of Lemmas 14.11 14.21 and 14.31 

Proof of Lemma\4-1[ First we show that for any map of the form u(x) := Aa + /(x) 



the 



incompressibility equation reduces to an ODE of the form f n 1 (r, C)§£( r > C) 



„n-l 



. In order to 



46 



see this, consider a local parameterization (si, . . . , s n -i) t- > C( s i> ■ • • > *n-l) °f 1 an d introduce 
polar coordinates of the form 

x = x(r, si, . . . ,s n -i) = a + rC(si, . . . ,s re _i), r > 0, (s x , . . . ,s„,_i) EflC E" -1 , 

D being some parameter space. The claim follows by observing that 



du du du _ . , 

— A — A • • • A = det Du x 

or os\ c*s n _i 



cbt <9x 



5x 



9r ds i dsn-i 



n-l 



det Du(x) C A /\ 



k=l 



ds k 



and 



du du du df . " A 1 / df . , dC\ .„_i 9/ 
A A • • • A = — C A A — C + /— = f — 

dr ds\ ds n -i dr \<9sfc ctefc / 5r 



n-l 



ca a 



fc=i 



From the above we find that u(x) := Aa + /(x)£ is incompressible provided f(r,£) n = r n + 
A(C) n , for some A : S™ -1 -> R. The definition in (j4"7T) . namely, f n = r n + |^</\ is obtained by 
imposing the boundary condition u(x) = Ax on dVt. Differentiating (|4.ip with respect to £ yields 



f n -\r, £)D c f(r, C) = r^q n - L (C)Dq(C, D c f(r, Q, Dq(Q : T c (! 



;n— In 



T^(S n ) being the tangent plane to § n at £. Identifying, in the usual manner, 



,n-l 



(0 



|0| /"-i(r,C) 

with a vector in (0 C K n , from /(x) = /(r(x), C( x ))> r ( x ) = l x — a l) C( x ) — | x _ a | 



(4.19) 



5/ 



r + (Dt) T D c f= d J-t + ^l, 



df 



dr 



+ 



i> q 71 ' 1 \Dq\ 



1 we obtain 

2 



n-l 



(4.20) 



with Dr = C and £>C = Since 1)11 = C®Df+fDC and {D£)-{C®Df) = {-((DQDf) = 0, 



using that |L>C| 2 = ^ and = £^ < 1 we find 



Du| 2 = | J D/| 2 + / 2 | J DC| 2 = (n-l)4 + 



9/ 



9r 



+ 



< fn- + l + 



(4.21) 



f 2 

The leading order term (vi + ^2)! log e| in the energy estimates will come from (n — 1)^, hence 



we need to write 



Du 



/n^l 



as plus a remainder (for which we do not require an exact expression, 



only an upper bound). To this end we bound a n — b n , with a 



Du 



/n-l 



and b = J ^ + by 



Du 



n — 1 r z 



< (a-6)|a n - 1 + --- + 6 n - i | <n 



|a 2 -6 2 | 



a + 6 



max{a, 6} 



n-l 



< 



n 



\iPcfVl 

n — 1 a + 6 



■max{a,6} n_1 < 



, n |(D c /)/r| 



n — 1 a 



max{a, 6} 



n-l 
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Prom f n = r n + j^q n and (j4TT9T) we find that 



f n v Q 

— = 1 + 



yn 



\n\ r»' 



/ > — T q, and \D c f\ < —r\Dq\. 



(4.22) 



As a consequence of (g^Hl , V^Ta > J^i, hence 

i , Z 2 ^ 



|Pc/)A| 2 < l^c/l 



and 



n — 1 r z ' 



(4.23) 



(we have used (|4.22|) to bound and (|4.2ip to bound max{a, 6}). Proceeding analogously, 

writing c = £ > 1 and b n — c n < n b2 b ~f 6 n ~ 1 < n(b 2 — c 2 )b n ~ l , we obtain 



1 , f 



v q 



+C 



+ 



1 + r 2 ) ~ 1+ |Q| r n ■ ' ~ V(n- l)f ' r" 



1 



< 



r n 



+ C 1 + 



i \ w— i 



. (4.24) 



Writing a n = 6" + (a n - 6 n ), equations g2D, (fOHl) . and <|Q2]> yield 
Du n 



„n-l 



< r 1 



n-l 



w\-- +c 1 



i \ n— 1 



i ; n 



+ C r |D?|r n - 2 1 + 



w 2 ^ 1 max{g n-1 , iDg]"- 1 } 

, „ , n — 1 

\n\— 



r n— 1 



l \ n—1 



< C r + 



L ^ n ._ 2 /g" max{g^g^pg| \ i> 
H r\Dq\r + tttt + -■ 

|n» vl n l N /r 



To finish the proof substitute both r and -^rlDol in -^r-\Dq\r n 2 with r + lul » . 



□ 



Proo/ o/ Lemma \JJ% Write C = cosfle + sin6>£, 9 G (0,7r), £ G S n_1 n (e)- 1 . By virtue of |(a + 
Q(C)C)-a| 2 = p 2 , 



g 2 + 2gC-(a-a) = p 2 -d 2 , q(6,£) = -dcos 9 + yj{p 2 - d 2 ) + d 2 cos 2 9. (4.25) 
Extending q to W 1 by g(x) = g(£(x)), C( x ) := rap an d differentiating with respect to x, we obtain 

(2q + 2C • (a - a)) P>g = -2g(P>C) T (a - a) = -2q^—^^-(de). 

|x| 

Our aim is to obtain bounds for q : S 71 " 1 — >• S n_1 and Dq(£) G (T^S n_1 )*. We can identify Dq 
with a vector in (C)" 1 " in the usual manner. From the relation Dq ■ e = Dq ■ ((D£)e) we know that 
Dq _L C and |Dg(x)| = |P>g( x )l for all x G S"" 1 . Thus, using that g 2 + 2<?C • (a - a) = p 2 - d 2 , 



\Dg(e,s)\ 



-2dq 2 (l -C®C)e 



g 2 + (g2 + 2gC-(a-S)) 



< 



2dq 2 sin 9 



< 



2dq 2 sin 9 



max{g 2 , (p — d)(p + d)} max{q 2 , p(p — d)} ' 
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this yields the bounds for \Dq\ in [T]). The fact that \q(C)\ < 2p for all C G § n 1 follows from 
g(C) = dist(a + g(C)C,a) < diam J3(a, p). Part |n|) is proved directly from the second equation in 
(|4,25p . considering that \/a + b < yfa + y/b, that p(p — d) < p 2 — d 2 < 2p(p — d) and that ^7 > 7 
for all 7 S (0, 1). Indeed, if cos# < then 

2d\ cos 6\ + y/2p{p - d) > d\ cos 6\ + yjp 2 - d 2 + Vd 2 cos 2 6 > q{9, £) 



> d\ cos 6\ + \J p 2 - d 2 > d\ cos 6\ + \fpjp - d) > p\ cos 0\ ^- + y 1 - - J . 
To prove [m]) , suppose that C " e = cos # > and rewrite (|4.25|) as 



□ 




Figure 10: Cone generated by S and a + pe (Lemma 14.3ft 



Proof of Lemma \4-3\ Call a := a + (p — d)e. Consider the (n — 2)-sphere S := {x € dB(a,p) : 
(x — a) • e = 0}. It is clear that contains the cone generated by a + pe (the 'right-most' point 
on dB(a.,p)) and S. Since the radius of 5 (the 'height') is given by h = y/d(2p - d) (see Figure 
[TO]) and the base measures d, the volume of the cone is a constant times d/i n_1 = d~^~ (2p — d)^~ . 
The value of the constant is obtained from 

101 > ^ n " 2(s "" 2) r ( p -^y/p 2 -(p-d) 2 ) n ~ l dxi = ^v^^r 1 -. 

□ 



4.4 Numerical computations 

The deformations depicted in Figure [6] are obtained by the alternative method of Dacorogna- 
Moser (constructive in nature and easier to implement, \20\ Sect. 4]). Following the notation in 

Theorem [3] (and restricting now to the case n = 2), let p{9) := \ R 2 + (v\ + v 2) WT^jn 1 1 where q{9) 
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denotes the parametrization of d(fli U 2 ) using polar coordinates, taking a* to be the origin. Let 
also < Ri < R 2 < R3 be such that B(a*,Ri) D fij U $7 2 and 7ri?§ = + t> 2 + nR 2 - Given 
parametrizations w(s,i) and v(s,t), (s,t) £ D := [1, v2] X [0, 27r] of {x : i?! < |x — a*| < R 2 } 
and of {y = Aa* + re %e : p{6) < r < -R3}, respectively, the strategy is to find an incompressible 
homeomorphism u : w(Q) — > v(Q) of the form 

u = vo0 2 o0 lO w, with </>i(s,t) = (h(s,t),t), <f> 2 (s,t) = (s,t + r}{s)/3(t)). 

Here ?] : [1, y/2] — > R is any function satisfying 

^ ??(s)ds = l, 77(0) = 77(1) = 0, 0<r/<l + e, j- |1 - 77(a)! ds <e 

for some e < mm j 2 ™"/ g , | , where /(s,t) = detDw(s,t) and #(s,i) = det Dv(s, i). The 
functions /3 and h are found by defining <7i(si,ii) := g(4> 2 (si,ti)) det Dcf> 2 (si, ti) and solving 

rV2 rt+^{a)l3(t) rV2 rt rh(s,t) rs 

/ / <7(a,r)drda= / f(s,t)dids, / gi ( Sl ,t)d Sl = f(s,t)ds 

Jl Jo Ji Jo Jl Jl 

for every fixed t € [0, 27r]. The solution is unique, and for v and w as in (|4.12p . it is such that 
fn <| _ *\<r l^ u l 2 — where C is an expression that might possibly go to infinity only if the 

tjge. do.nail is too narrow, m ore precis*, if ,V » - 1, (recall that 

2 

l^yo^l — 1 is of the order of 1 — S, equations (|4.1U|) and (|4.15|) ). In our computations we choose 

,2 
m; 

|J2iun 2 | 

5 Proof of the convergence result, Theorem |4] 

We follow the strategy of Struwe [77J to prove that sup e ||u e ||wi,p(n e ) < oo for all p < n. Fix e > 0, 
call £>o := \JiLi B £ (a.i j£ ), to := r(Bo) = m^-, and let : t > to} be the family obtained by 

applying Proposition 13.21 to Bq. Define p = sup{t > to : \jB(t) C SI} and write Ck ■= \jB(rk) \ 
U B(rfe+i), rfc := 2~ k p. By using Holder's inequality, then comparing the lower bound of Proposition 
13.41 to the upper bound, we find that for every p < n 



Ri = <7ma X = 2pi - d5 and i? 2 such that tt(R% - R\) = 2{vi + v 2 ) ( r^fWi ~ 1 ) • 



/ |L>u £ | p dx< C(n,p)p n - p 2-( n - p ) fc - / ; Ue dx-Vt; le lo 
Jc k \nJn E y/n-1 

/ m 

< cp n -n~ [n -^ k ( |n| + 



— 



to 
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Adding over k we find that 



|-Du e | p dx < Cp 



n-p 



<C(p 



n—p 




It can be seen (as in the proof of Proposition II. ip that p > \ dist({ai j£ , . . . , a mi£ }, dQ). Hence, in 
order to prove that sup £ ||Z)u e ||iP < oo it only remains to show that YliLi v i,e ^ s uniformly bounded. 
Choose r > e such that the balls £>(aj j£ , r) are disjoint and r E R &i e for alH = 1, . . . , m. By Propo- 
sition [2TT| the topological images E{sti j£ , r; u £ ) are disjoint, contained in -6(0, Hu^^oo^)) (because 
E(ai >e , r; u E ) is the region enclosed by u(dB(a.i j£ , r))), and such that D(aj )£ , e; u e ) C E{&i i£ , r; u £ ). 
Therefore 



in 

£ 

i=l 



£(ai !£ ,£;u e )| < 



[J £(a ii£ ,r;u e ; 



< 0J r , 



U F 



Since we are assuming that sup £ ||u e ||x,cx>(Q e ) < oo, we obtain that sup £ ||u e || w -i, P (Q e ) < oo, as desired. 



For the existence of a limit map and for the convergence in Wj '"(fi \ {ai, . . . a m }, W 1 ), let 5 > 
be small, assume that |aj e — a^| < 5/2 for alH = 1, . . . , m, and consider the following energy bound, 
obtained again by comparing (|1.16p with the lower bound of Proposition 13.41 (applied to s = 5/2) 



1 

n Jn\[JB(8/2) 



Du 



n 



1 



dram SI / x - 
0/2m V f-f 

1 = 1 ' \ 2=1 



Since r{B{8/2)) = 5/2, it follows that {u e J jeN is bounded in W 1>n (n \ U=i ^(aj), M. n ). From 
this, and since 5 > is arbitrary, the existence of u and of a convergent subsequence follows by 
standard arguments (see, e.g., [72] or [37J): inductively take succesive subsequences of {u £j }j 6 N 
(for some sequence 5k — > 0) converging weakly in W 1,n (Q \ \J^ =1 Bs k (ai),M. n ). Choose then a 
diagonal sequence {u £fc }fc 6 ^ converging weakly in W 1,n (Q. \ {J//L x B$(ai), K n ) for every 5 > 0, to 
some u e W^(n \ {a l5 . . . a m }, E«). 



Since sup £ ||u e 



< oo for all p < n, the maps u £ can be extended, by multiplying them by 



suitable cut-off functions ip £ , inside the holes B(a.i y£ , e), in such a way that sup £ ||'i/; e u E ||vK 1 .p(n) < °°- 
It is easy to see that any weakly convergent subsequence of {ip £k u £k }k^N must converge to the limit 
map u defined above; this proves that u G W 1,P {Q., R n ) for all p < n. 

By the classical result of Reshetnyak |61|, Thm. 4] and Ball Cor. 6.2.2], cof Du £k — ^ cof Du 

n 

in L 1 , ;; 1 (Sl\{ai,...,a m },]R nxn ). By the definition of DetDu in ([53]) . and since {DetDu e } e>0 is 
bounded as a sequence in the space of measures (Det Du £ = £ n L £ , by hypothesis), it follows 
that Det Du coincides with C n in \ {ai , . . . , a. m }, and that Det Du £ — Det Du in f2 \ {ai , . . . , a m } 
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in the sense of measures. Moreover, by [70j Lemma 3.2] (applied to Q \ U^Li B(ai, 5) instead of $7), 
we obtain that detDu(x) = 1 for a.e. x G Q \ {ax, . . . ,a m }. 

From Definition [5] and from the proof of |37^ Lemma 4.2] it follows that the limit map u satisfies 
condition INV. Proposition 12.21 then implies that Det Du = C n + Y2i=i c i^a.i f° r some coefficients 
Cj G R, and the proof of the same proposition also shows that 

= u n r n + ^ v j,e 



- [ u £ • (cof Du £ )v dH^ 1 = u n i 

-f u (cofDuWd?^ 1 = u} n r n + V c 7 - 
n JdBU„ r ) . ~ s 



j:a jiS &B(ai,r) 



j:a J eB(a 1 ,r) 



for a.e. r > such that dB{s^, r) C O (note that if aj = a.j for some i ^ j, then the choice of the 
coefficients c% is not unique). By standard arguments, for every 5 > there exists r < 5 such that 

re 

u £k — > u uniformly on <9£?(a.j,r) and cof Du E| . — cof Du in L n - 1 (dB(&i,r)) (passing, if necessary, 
to a subsequence that may depend on r). Taking, first, the limit as e — > 0, then the limit as r — > 0, 
we obtain that Det Du = C n + YliLi u i<W 

Consider now the case of two cavities. Set a e := ai ' E ^ a2,e , d e := |a2 )E — ai j£ |. 

i) Suppose that v% > i>2 > and d = \a2 — ai| > 0. By Lemma 13.51 we have that for all r > e 

n n re 

|£(a i!£ ,r;u e )|£(£K e ,r;u e ))^ - ^(a^, e; u e )|D(E(a ii£ , e; u £ ))^ < 2^£±±uv n , 
hence, by (|3.5p . for all a G (0, 1) and all R < min{|, dist({ai, a2}, dtt)} we have that 



Jr2 £ n 



Du(x) 



loge| 



> 



Yli=lif £ + J £ Q J IdB(a i>£ ,r) n 



Du(x) 



dTT^dr 



loge| 



> E(^ 1 °ff-f ) - f 



loge| 



Combining this with (|1.16p we obtain 



X — A re 

^^(^(a^e.-u,))"- 1 < 



\n\+v ljj + v 2 , j )(C 2 + \og^) 



i=l 



Ci | log e 



l-a\ 



+ Ce ari . 



Therefore, as e — > 0, Z?(£ , (aj i£ , e; u e )) — > (i.e., u e tends to create spherical cavities). 

As mentioned before, for every 5 > there exists r < 5 such that u e |g^( a . r ) converges uni- 
formly, for each i = 1,2, to u|gs(a l ,r) (passing to a subsequence, if necessary). By continuity 
of the degree, this implies that iniT(u, a^) is contained in E(ai,r; u £ ) for sufficienty small e. In 
particular, by definition of Vi £ and Proposition 12.21 



|£ , (ai,r;u £ )Aim T (u,a i )| = \E(ai,r;u e )\ - | im T (u,aj)| = (vi >e + io n r n ) - Vi 
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On the other hand, B(ai e ,e) C B(a.i, r) for sufficiently small e. By Proposition 12 . 1 1 this implies 
that £ , (aj e , e; u e ) C E(ai,r;u £ ), so, proceeding as in the proof of Proposition 12.21 we obtain 

|£(a ii£ ,e;u e )A£(ai,r;u £ )| = T)etDu(B(a i ,r)\B(a iiE ,e)) = \B(a il r)\B(a iill ,e)\ < uj n 5 n . 

Thus, 

limsup \E(a it£ , e; u e ) A im T (u, aj)| (5.1) 

< limsup(|-E(a ii£ ,e;u e )A£;(a i ,r;u £ )| + |£(aj,r;u £ )Aim T (u,aj)|) < 2uj n 5 n 

for all 5 > 0, that is, the cavities formed by u £ converge to the cavities formed by u. 
It remains to prove the estimate for |a2 — ai| in terms of |f2| , diamil and the cavity volumes, 
assuming that V! + v 2 < 2 n w n (dist( ai ± ai ,<90)) r \ Let R > be such that vi >s + v 2 , £ < oj n (2R) n 
and B(a £ ,R) C SI for every sufficiently small e. Suppose first that 

^ <U^-Y\ (5.2) 



V\ + V 2 2 n \Vi + V 2/ 

Since v ^ v < 1, this implies, in particular, that «i )£ + t> 2j£ > co n (2d £ ) n for every small e. As 

i j_ 

a consequence, ^ > 1 and ( ^n^d™ ) < (<^) ^ * na ^ * ne m ffii m u m at the end of 

Theorem [T] is attained at ( ~2^uf~3* ) cann °t be attained at ^- since d £ — >■ (i > 0). By 

Theorem [T] and (11 . 16f) . 



^2,e \ "- 1 ^n^ £ 



C,l I " ) _ — , log ^±^<Mb ^ dx- Ke + „ 2 ,.)log 
^ \Ul )£ + «2,e/ Vl,e + V2,eJ 2 n UJ n d £ l V\,e + V 2)£ 

_ / |n| , a;„(diamO) n \ _ / 101 , w n (diam$7) n 

< C 2 1 + U + log — — < C 2 1 + U + log 



Vl,s + V 2 , £ UJ n R n J V v l,s + ^2,e Ul,e + V 2>£ 

(in the last step we use that uj n R n > Vl ' £ + V2 ' e , by the choice of R). If O holds then the 
factor in front of the logarithm is positive for e > small; taking the limit we obtain that 
agL>2-»F(J2, with 

F(fi, := ex P (-C (l + -H— + log """" ) / f . (5.3) 
y V vi+v 2 Vi+V 2 )l \v 1 +v 2 J J 

If (|5.2p does not hold, we still have that > CF(fi, «i, ?; 2 ) for some constant C(n). To see 

this, recall that v\ + v 2 < 2 n ui n dist( ai ' 2 " a2 , <9f2)) n < w n (2 diam 0) n (by hypothesis), hence 

F(0, Vl , U2 )<exp( -C(l + n|log2|) / ( V"' ) < f ~ Z~ 1 ^ / C(l + n|log2|) 

(we have used that > - for all cc > 0). The proof is completed since the above implies that 

n 

^ > 2"» =► ^^>2-C7(l + n|log2|)F(0,, 1 ,, 2 ). 



«1 + V 2 \Vl+V 2 J V 1 + ?J 2 
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ii) Suppose that v\ > v 2 = 0. Applying Proposition 13.21 to Bo := {B e (ai >e ), £? £ (a 2j£ )} we obtain 
B(t) = {J3(a 1>e ,V2),J3(a 2je ,</2)} for t G (2e,d £ ), and B(t) = {S(a e ,t)} for t > d £ . We 
claim that if R < | dist({ai i£ , a 2;£ }, dQ,) then lj£?(_R) C O. Indeed, if it! < d £ , this holds 
automatically. If R > d s , then 

^ < dist(a 1)E , dfl) < y + dist(a e , 917) < ^ + dist(a £ , 90) => £(a £ , ft) c O. 

Therefore, by Proposition 13.41 and Lemma 13.51 for every a G (0, 1) 

|-E(ai, e ,e; u £ )|£>(£(ai i£ ,e; u E )) log — 



< 



/ 1 




In, n 


y/n-1 



dx - (v 1)B + V2, e ) log |- + 2^T^±i( U2 , £ + Wn e a ") log 



By virtue of (|1.16p and again Lemma 13.51 

^(im^u, ai)) - 1 < 2^1 2±i lim(^ 2i£ + w n e an + |£(a lj£ , e; Uj -)A im T (u, ai)|). 
Proceeding as in (15. ip we find that 

limsup|£(ai i£ ,e;u £ )Aim T (u,ai)| < 2(v 2 + uo n r n ) 

for arbitrarily small values of r > 0, proving that imx(u, ai) is a ball. 

Suppose that v\ > -u 2 > and ai = a 2 . Let R > be such that B(a £ ,R) C f2 for all j G N. 
Since limd £ = |a 2 — ai| = 0, (|3.6p and f 1 1 . 1 6 1) imply that 



f d R \E(a e ,r; u £ )\D(E(a £ ,r; u £ )) ^ + t*) (l + log ^ 

bmsup— * — < C — - — = 0. 

£ ^ log d £ lim log d £ 

£->0 

For a G (0, 1) fixed and e small B(a e , d £ ) C B(a e , d £ ) C O. By Lemma [331 for all r G (d £ , d £ 



|£(a e ,r; u £ )|Z)(£(a £ , r; u £ )) - |£(a £ , d e ; u £ )|L>(£(a £! d £ ; u £ )) » 
Dividing \E(a £ ,d £ ;u £ )\D(E(a £ ,d £ ;u £ ))^^ by log^" 1 we obtain 



n 71+1 

<2— _L^ Wn d^. 



n — 1 



ti n 72 —I— 1 

limsup|£(a £ ,d £ ;u £ )|Z)(£(a £ ,4;u £ ))"-i < lim sup 2 ~ -uj n df n = 0. (5.4) 

Proceeding as in (15. ip . it can be proved that 

limsup I im T (u,ai)A£ , (a £ ,d £ ;u £ )| < limsup(t>i j£ + v 2 , E ) - | im T (u,ai)|. (5.5) 

£^0 £^0 

Because of the continuity of the distributional determinant, | imT (u, ai ) | = V\ + u 2 , hence 
D(im T (u,ai)) = (by Lemma E30), and 
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In order to prove that at least one of the limit cavities must be distorted, we proceed as in 

the proof of TheoremQ]by applying Proposition 11.21 to E\ = E(a% E , e; u e ), E 2 = E(a 2e ,e;u e ), 

x x 

and E = E(a e ,d e ;u e ). Again we define /3 2 ) := (/3f + /3 2 ") n — {/3\ + /3 2 ) an d note that it 
is increasing in its two variables. It is easy to see that 

(\Et\^ + \E 2 \^) n -\E\ UJ n d n e e^o 

> 1 ; r > 1. 



+ \E 2 \") n - \E 1 U E 2 \ 9(vi,e,V 2 ,e) 

Therefore, 

,. . £ \E\D{E)^ + \Ei\D{E 1 )^ + \E 2 \D{E 2 )^ ^ ^ f v 2 
hm mi — ; — — > C 



e^o \E\ + \Ei U E 2 \ \vi+v 2/ 

Property (|1.17j) follows from (|5.4p . On the other hand, (|3.6p . (|1.16|) . and Lemma 13.51 imply 
that 



f /■min{-^,e"} / n jn \ j 



for every fixed a £ (0, 1). Hence, 



£->0 



. , . diamfi . . . 

< [Vl,e + V 2tE ) log — + 6 + U 2 , e + |0| 



diamSl 
C log — + 1 + 



, 4 , a-l\\ V A/2 Vl+V 2 , 

hmsup mm ^ log — ,loge > < 7 s s ^ 

2e J / Um inf f vi,eD(E 1 )^+v 2 , £ D(E 2 )^ _ ^ 



By virtue of (|1.17p . and since | loge| — > oo, we conclude that limsupci e /e is finite. 

e-s>0 
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